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GRANT#:  N00014-05-1-0009 

PRINCIPAL  INVESTIGATORS:  Diola  Bagayoko,  Ph.D.,  Southern  University  System 
Distinguished  Professor  of  Physics,  Project  Director,  and  G.  L.  Zhao,  Ph.D.,  Associate  Professor 
of  Physics,  SUBR,  Co-Principal  Investigator 
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Ballston  Centre  tower  One,  800  North  Quincy  Street,  Arlington,  VA  22217-5660 

INSTITUTION :  Southern  University  and  A  &  M  College  in  Baton  Rouge  (SUBR) 

GRANT  TITLE:  Predictive  Computations  of  Properties  of  Wide-Gap  and  Nano- 
Semiconductors 

AWARD  PERIOD:  October  1,  2004  to  September  30,  2005 

OBJECTIVE:  As  per  the  initial  proposal,  the  aim  of  this  project  is  (a)  to  produce  new  scientific 
knowledge  on  wide-gap  and  nano-semiconductors,  with  specific  materials  including  GaN,  InN, 
AIN,  and  carbon  nanotubes,  (b)  to  develop  analytical  and  computational  techniques  and 
computer  codes  for  wider  applications  in  descriptive  and  predictive  investigations  of  the  above 
materials  and  others,  and  (c)  to  provide  research  training  for  undergraduate  and  graduate 
students,  including  minority  students  at  SUBR. 

APPROACH:  Our  technical  approach  or  method,  as  per  the  proposal,  consists  of  the  following. 
Like  most  previous  calculations  of  electronic  and  related  properties  of  semiconductors,  we 
employ  a  density  functional  theory  (DFT)  potential.  This  potential  is  often  selected,  by  design, 
from  the  local  density  approximation  (LDA)  to  the  general  DFT.  We  also  utilize,  like  previous 
works,  the  linear  combination  of  atomic  orbital  (LCAO)  formalism  in  a  Rayleigh-Ritz  variational 
procedure  to  perform  self-consistent  calculations.  Our  calculations,  to  date,  are  for  zero 
temperature  and  are  non-relativistic. 

The  key  distinction  between  our  method  and  those  of  all  previous  calculations  by  other 
groups  resides  in  our  adherence  to  the  Bagayoko,  Zhao,  and  Williams  (BZW)  method  in  the 
implementation  of  the  LCAO  formalism.  References  clearly  identified  in  the  attached  articles 
fully  describe  this  method.  Essentially,  this  method  resolved  the  long-standing  underestimation, 
by  theory,  of  measured  energy  or  band  gaps  of  finite  (atoms,  clusters,  molecules)  and  infinite 
(semiconductors  insulators)  systems.  As  explained  in  the  attached  articles  and  references  therein, 
BZW  first  identified  a  basis  set  and  variational  effect  that  is  responsible  for  an  unphysical 
lowering  of  the  low-lying,  unoccupied  energy  levels  or  bands  of  materials.  This  unphysical 
lowering  is  shown  to  be  a  mathematical  artifact  stemming  from  the  Rayleigh  theorem.  BZW 
subsequently  introduced  a  rigorous  method  of  performing  variational  calculations  in  such  a  way 
that  the  above  unphysical  lowering  of  unoccupied  levels  or  bands  is  totally  avoided.  Let  us 
underscore  that  the  underestimation  of  energy  and  band  gaps  is  a  trivial  consequence  of  the 
above  unphysical  lowering  of  unoccupied  levels  of  bands. 
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The  BZW  method  generally  entails  the  performance  of  three  (3)  or  more  self-consistent 
calculations.  The  first  calculation  employs  the  minimal  basis  set  for  the  system  and  progressively 
augment  it  in  subsequent  calculations.  The  occupied  energy  levels  or  bands  of  a  calculation  the 
next,  consecutive  calculation  whose  basis  set  is  larger.  As  explained  in  detail  in  the  articles, 
successive  calculations  have  different  basis  sets:  the  basis  set  for  the  latter  is  obtained  by 
augmenting  that  of  the  former  with  the  next  excited  state  orbital  an  atomic  or  ionic  species  that  is 
present  in  the  system  under  study  (i.e.,  In  and  N  for  InN).  This  process  continues  until  the 
occupied  energies  are  found  to  be  identical  for  two  consecutive  calculation  N  and  (N+l).  Then, 
the  results  of  calculation  N  are  the  physical  ones  for  the  system  under  study  (and  not  those  of 
calculation  N+l).  This  criterion  for  stopping  the  process  rests  trivially  on  the  facts  that  (a) 
density  functional  theory  is  a  ground  state  theory  and  (b)  only  the  wave  functions  of  the  occupied 
states  are  including  in  the  construction  and  reconstruction  of  the  charge  density  (and  hence  the 
potential  and  the  Hamiltonian)  in  the  self-consistency  process. 

Once  the  correct  energy  levels  or  bands  are  obtained,  then  we  proceed  to  calculate  a  host 
of  properties  that  include  the  energy  or  band  gaps,  optical  transition  energies,  the  total  energy, 
the  total  and  partial  densities  of  states,  the  effective  masses  of  charge  carriers  (for  application  in 
transport  studies),  and  the  optical  properties  (i.e.,  the  complex  dielectric  functions).  Of  course, 
the  total  energy  curves  are  utilized  to  compute  the  bulk  modulus. 

In  addition  to  the  above  technical  approach,  we  should  note  that  our  programmatic 
implementation  also  adhered  to  the  proposal.  In  particular,  we  have  hired  Dr.  Hua  Jin  as  the 
Postdoctoral  Associate.  She  has  been  performing  excellently,  not  only  in  carrying  out  research 
tasks,  but  also  in  assisting  the  investigators  in  their  training  of  students  on  the  complex  software 
package  utilized  in  this  work.  A  graduate  and  undergraduate  students,  mostly  African  American, 
have  also  been  employed.  SUBR  is  a  Historically  Black  College  and  University  (HBCU).  The 
effort  levels  of  the  investigators  (D.  Bagayoko  and  G.  L.  Zhao)  are  as  per  the  proposal. 

ACCOMPLISHMENTS:  As  summarized  in  the  abstract  of  this  report,  we  have  published  six 
(6)  articles,  made  seven  (7)  technical  presentations,  two  of  which  were  at  an  international 
conference.  A  seventh  article  is  pending.  To  reduce  redundancy,  we  refer  the  reader  to  the 
following  full  listing  of  the  articles,  with  referencing  information.  We  also  list  the  presentations 
below.  The  subtitles  below,  for  accomplishments,  are  keyed  to  the  objectives  enunciated  above. 

Publications 

“Structural,  Elastic,  and  Electronic  Properties  of  Deformed  Carbon  Nanotubes  under  Uniaxial 
Strain.  ”  A.  Pullen,  G.  L.  Zhao,  D.  Bagayoko,  and  L.  Yang.  Physical  Review  B  71,  205410. 

“Ab-initio  Simulations  of  the  Growth  and  Structural  Properties  of  Short  Carbon  Nanobells.”  G. 
L.  Zhao,  D.  Bagayoko,  and  E.  G.  Wang.  Accepted  for  publication  in  the  Proceedings  of  the  2005 
China  Conference  on  Nanoscience  and  Technology,  Beijing,  China,  July  2005. 

“A  Universal  Relation  Between  the  Densities  of  States  Near  van  Hove  Singularities  and  the 
Effective  Electron  Masses  in  1-Dimensional  Semiconductors.  ”  G.  L.  Zhao  and  D.  Bagayoko. 
Accepted  for  publication  in  the  Proceedings  of  the  2005  China  Conference  on  Nano-Science  and 
Technology  ,  Beijing,  China,  July  2005. 
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“LDA  and  LCAO-BZW  Description  of  Electronic  Properties  ofWurtzite  Zinc  Oxide  (w-ZnO).  ” 
Diola  Bagayoko,  Lashounda  Franklin,  and  G.  L.  Zhao,  accepted  for  publications  in  Proceedings 
of  the  2005  National  Conference  of  the  National  Society  of  Black  Physicists  (www.nsbp.org). 
Orlando,  Florida. 

“Predictions  of  Electronic,  Structural,  and  Elastic  Properties  of  Cubic  InN.  ”  Diola  Bagayoko, 
Lashounda  Franklin,  and  G.  L.  Zhao,  Journal  of  Applied  Physics  96,  4297-4301,  2004. 

“ Density  Functional  Band  Gap  of  Wurtzite  InN.  ”  Diola  Bagayoko  and  Lashounda  Franklin, 
Journal  of  Applied  Physics,  97,  123708,  2005. 

“Re-examination  of  the  Ab-initio  Calculation  of  the  Electronic  Structures  of  ZnSe,  Ge,  and 
GaAs.  ”  G.  L.  Zhao,  L.  Franklin,  and  D.  Bagayoko.  Submitted  to  Physical  Review  B,  2005. 


Presentations 

June  9-11,  2005,  Beijing,  China.  China  International  Conference  on  Nanoscience  and 

Technology.  “Ab-initio  Simulations  of  the  Growth  and  Structural  Properties  of  Short  Carbon 
Nanobells.  ”  G.  L.  Zhao,  D.  Bagayoko,  and  E.  G.  Wang. 

June  9-11,  2005,  Beijing,  China.  China  International  Conference  on  Nanoscience  and 

Technology.  “A  Universal  Relation  Between  the  Densities  of  States  Near  van  Hove  Singularities 
and  the  Effective  Electron  Masses  in  1 -Dimensional  Semiconductors.”  G.  L.  Zhao  and  D. 
Bagayoko. 

March  25,  2005.  Los  Angeles,  CA.  March  Meeting  of  the  American  Physical  Society  (APS). 
“Structural,  Elastic,  and  Electronic  Properties  of  Deformed  Carbon  Nanotubes  under  Uniaxial 
Strain."  A.  Pullen  SUBR  &  (Caltech),  G.  L.  Zhao,  D.  Bagayoko,  and  L.  Yang  (NASA),  Bull. 
APS,  Vol.  50,  No.  1,  Page  1420  (2005). 

March  23,  2005.  Los  Angeles,  CA.  March  Meeting  of  the  American  Physical  Society  (APS). 
“Re-examination  of  Ab-initio  Calculation  of  the  Electronic  Structure  of  Zn  Se,  Ge,  and  GaAs.” 
G.  L.  Zhao,  L.  Franklin,  and  D.  Bagayoko,  Bull.  APS,  Vol.  50,  No.  1,  Page  1073  (2005). 

March  22,  2005.  Los  Angeles,  CA.  March  Meeting  of  the  American  Physical  Society  (APS). 
“True  LDA  Band  Gaps  ofWurtzite  and  Cubic  Indium  Nitride  (w-InN  and  c-InN).  ”  D.  Bagayoko, 
G.  L.  Zhao,  and  L.  Franklin.  Bull.  APS,  Vol.  50,  No.  1,  Page  617  (2005).  Audience: 
Approximately  70  physical  science  researchers  and  graduate  students. 

February  18,  2005.  Orlando,  Florida,  Disneyland.  2005  National  Conference  of  the  National 
Society  of  Black  Physicists  (NSBP)  and  of  the  National  Society  of  Hispanic  Physicists  (NSHP). 
“Local  Density  Functional  Description  of  Electronic  Properties  ofWurtzite  Zinc  Oxide  (ZnO).” 
D.  Bagayoko,  G.  L.  Zhao,  and  L.  Franklin.  Audience:  17  faculty  members,  graduate  students, 
and  federal  lab  researchers. 
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February  17,  2005.  Orlando,  Florida,  Disneyland.  2005  National  Conference  of  the  National 
Society  of  Black  Physicists  (NSBP)  and  of  the  National  Society  of  Hispanic  Physicists  (NSHP). 
“A  Competitive  Edge  for  Recruitment:  The  Versatility  and  Wonders  of  Physics.  ”  D.  Bagayoko. 
Audience:  20  faculty  members,  graduate  students,  and  federal  lab  researchers. 

Code  Development 

Another  major  component  of  our  accomplishments  in  this  reporting  period  stems  from  the 
development  and  refinement  of  the  computer  program  package  for  wider  utilization.  We  have 
developed  a  precise  and  concise  user-manual  that  we  are  currently  field-testing  for  effectiveness. 

Research  Training  of  Students 

We  have  trained  four  (4)  minority  undergraduate  students  in  the  performance  of  the  sophisticated 
calculations  this  project  entails.  These  students  are  physics  majors  at  SUBR.  One  of  the  students 
is  the  first  author  on  the  Physical  Review  article  on  electronic,  elastic,  and  related  properties  of 
single  walled  carbon  nanotubes. 

CONCLUSIONS:  We  have  implemented  in  a  successful  manner,  as  gauged  by  the  publications 
and  technical  presentations,  this  project  focused  on  “predictive  calculations  of  properties  of 
wide-gap  and  nano-semiconductors.  Materials  studied  to  date  include  GaN,  ZnO,  wurtzite  InN, 
cubic  InN,  and  several  carbon  nanotubes.  Our  results,  in  contrast  to  previous  DFT  or  LDA 
calculations  by  other  groups,  are  in  excellent  agreement  with  experiment.  The  significance  of 
these  results  is  underscored  below  and  in  the  abstract. 

SIGNIFICANCE:  Our  articles  on  electronic,  structural,  optical,  and  other  properties  of 
semiconductors  and  carbon  nanotubes  (Please  see  the  attached  list  of  publications)  have 
profound  implications,  even  for  the  study  of  nuclei.  These  papers  resolved  a  problem  dating  back 
to  the  beginning  of  quantum.  Indeed,  we  solved  the  woeful,  theoretical  underestimation  of  the 
band  gaps  of  semiconductors  and  insulators  as  compared  to  measured  values.  The  resolution  of 
this  problem  has  some  serious  implications  for  the  design  and  fabrication  of 
semiconductor-based  and  nano-devices.  A  collateral  benefit  stems  from  the  fact  that  the 
method  of  solution  also  applies  to  the  prediction  of  nuclear  energy  levels  in  the  shell  model, 
opening  the  way  to  theoretical  explorations  of  possibilities  of  a  population  inversion  and  the 
actual  construction  of  a  “gamma  ray  amplification  by  stimulated  emission  of  radiation  (graser)” 
device.  We  believe  that  no  comments  are  needed  on  commercial  and  other  applications  of  a 
graser. 

PATENT  INFORMATION:  Not  Applicable 
AWARD  INFORMATION:  Not  Applicable 

ATTACHMENTS:  FULL  TEXT  OF  PUBLICATIONS  LISTED  ABOVE 
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“Structural,  Elastic,  and  Electronic  Properties  of  Deformed  Carbon  Nanotubes  under  Uniaxial 
Strain.  ”  A.  Pullen,  G.  L.  Zhao,  D.  Bagayoko,  and  L.  Yang.  Physical  Review  B  71,  205410. 

“Ab-initio  Simulations  of  the  Growth  and  Structural  Properties  of  Short  Carbon  Nanobells.  ”  G. 
L.  Zhao,  D.  Bagayoko,  and  E.  G.  Wang.  Accepted  for  publication  in  the  Proceedings  of  the  2005 
China  Conference  on  Nanoscience  and  Technology,  Beijing,  China,  July  2005. 

“A  Universal  Relation  Between  the  Densities  of  States  Near  van  Hove  Singularities  and  the 
Effective  Electron  Masses  in  1 -Dimensional  Semiconductors.”  G.  L.  Zhao  and  D.  Bagayoko. 
Accepted  for  publication  in  the  Proceedings  of  the  2005  China  Conference  on  Nano-Science  and 
Technology  ,  Beijing,  China,  July  2005. 

“LDA  and  LCAO-BZW  Description  of  Electronic  Properties  ofWurtzite  Zinc  Oxide  (w-ZnO).” 
Diola  Bagayoko,  Lashounda  Franklin,  and  G.  L.  Zhao,  accepted  for  publications  in  Proceedings 
of  the  2005  National  Conference  of  the  National  Society  of  Black  Physicists  (www.nsbp.org). 
Orlando,  Florida. 

“Predictions  of  Electronic,  Structural,  and  Elastic  Properties  of  Cubic  InN.  ”  Diola  Bagayoko, 
Lashounda  Franklin,  and  G.  L.  Zhao,  Journal  of  Applied  Physics  96,  4297-4301,  2004. 

“ Density  Functional  Band  Gap  of  Wurtzite  InN.  ”  Diola  Bagayoko  and  Lashounda  Franklin, 
Journal  of  Applied  Physics,  97,  123708,  2005. 

“Re-examination  of  the  Ab-initio  Calculation  of  the  Electronic  Structures  of  ZnSe,  Ge,  and 
GaAs.  ”  G.  L.  Zhao,  L.  Franklin,  and  D.  Bagayoko.  Submitted  to  Physical  Review  B,  2005. 


PHYSICAL  REVIEW  B  71,  205410  (2005) 


Structural,  elastic,  and  electronic  properties  of  deformed  carbon  nanotubes  under  uniaxial  strain 


A.  Pullen,1  G.  L.  Zhao,1’*  D.  Bagayoko,1  and  L.  Yang2 
1 Department  of  Physics,  Southern  University  and  A  &  M  College,  Baton  Rouge,  Louisiana  70813,  USA 
2Eloret,  NASA  Ames  Research  Center,  MS230-3,  Moffett  Field,  California  94035,  USA 
(Received  7  May  2004;  revised  manuscript  received  20  January  2005;  published  25  May  2005) 


We  report  structural,  elastic,  and  electronic  properties  of  selected,  deformed,  single-wall  carbon  nanotubes 
under  uniaxial  strain.  We  utilized  a  generalized  gradient  approximation  potential  of  density  functional  theory 
and  the  linear  combination  of  atomic  orbital  formalism.  We  discuss  bond-lengths,  tubule  radii,  and  the  band 
gaps  as  functions  of  tension  and  compression  strain  for  carbon  nanotubes  (10,  0),  (8,  4),  and  (10,  10)  which 
have  chiral  angles  of  0,  19.1,  and  30  deg  relative  to  the  zigzag  direction.  We  also  calculated  the  Young’s 
modulus  and  the  in-plane  stiffness  for  each  of  these  three  nanotubes  as  representatives  of  zigzag,  chiral,  and 
armchair  nanotubes,  respectively.  We  found  that  these  carbon  nanotubes  have  unique  structural  properties 
consisting  of  a  strong  tendency  to  retain  their  tubule  radii  under  large  tension  and  compression  strains. 
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Single-wall  carbon  nanotubes  (SWCNTs)  can  be  viewed 
as  rolled-up  graphene  sheets  that  have  a  diameter  at  the  order 
of  1  nm.  They  have  properties  such  as  high  current  density, 
high  elasticity,  and  stiffness  unparalleled  by  other  materials. 
Their  potential  applications  range  from  that  in  building  sky¬ 
scrapers  and  elevator  cables  to  the  ones  in  very  tiny  electrical 
circuits  and  machines.1  These  materials,  however,  are  too 
small  for  many  conventional  measurements.  This  situation 
underscores  the  possible  importance  of  theoretical  studies, 
including  the  one  reported  here  that  focuses  on  structural, 
elastic,  and  electronic  properties  of  selected  single-wall  car¬ 
bon  nanotubes  under  uniaxial  strain. 

In  the  last  several  years,  tight-binding  calculations  have 
been  extensively  used  to  study  the  structural  and  electronic 
structure  of  carbon  nanotubes.2-6  Tight-binding  approxima¬ 
tions  based  on  the  symmetry  of  the  honeycomb  lattice  of 
graphite  predicted  that  SWCNTs  could  be  semiconducting  or 
metallic  depending  on  their  chirality  (n,m).  The  tight- 
binding  model  has  been  able  to  provide  good  estimates  of  the 
basic  electronic  structure  of  SWCNTs.  However,  curvature- 
related  effects  and  the  hybridization  of  different  electronic 
states  of  graphite  could  lead  to  structural  and  electronic  prop¬ 
erties  that  are  substantially  different  from  the  result  of  tight- 
binding  calculations.  Zigzag  (n,0)  (where  n  is  a  multiple  of 
3)  SWCNTs  which  were  predicted  to  be  metallic  from  tight- 
binding  calculations  were  found  to  have  small  energy  gaps.7 
Several  theoretical  groups  have  studied  the  elastic  properties 
of  carbon  nanotubes.  Their  approaches  include  simulations 
with  realistic  many-body  potentials,8  the  empirical  force- 
constant  method,9  tight-binding  formalisms,6,10  pseudopoten¬ 
tial  calculations  with  local  density  approximation 
potentials,11  and  Bom’s  perturbation  technique  within  a  lat¬ 
tice  dynamics  model.12 

The  aim  of  this  work  is  therefore  to  study  the  aforemen¬ 
tioned  properties  of  SWCNTs,  utilizing  ab  initio  quantum 
computations.  We  specifically  report  on  structural,  elastic, 
and  electronic  properties  of  SWCNTs  under  uniaxial  strain. 
Recent  tight-binding  calculations  have  led  to  values  of  bond 
lengths  and  radii,  band  gaps,  and  Fermi  levels  as  functions  of 
strain.6  The  dependence  of  tight-binding  results  upon  the  se¬ 


lected  parameters  partly  suggests  the  present  work  that  em¬ 
ploys  an  ab  initio  approach. 

Three  main  features  characterize  our  first-principles  com¬ 
putational  method.  The  first  of  these  features  is  the  choice  of 
the  potential.  We  utilized  the  generalized  gradient  approxi¬ 
mation  (GGA)  potential  of  Perdew  and  Wang.13  This  density 
functional14-16  potential  goes  beyond  the  local  density  ap¬ 
proximation  (LDA).  We  also  performed  LDA  calculations 
for  carbon  nanotube  (10,  0).  The  second  feature  of  our 
method  stems  from  employing  the  linear  combination  of 
atomic  orbital  (LCAO).  The  third  and  distinctive  feature  of 
our  work  resides  in  our  use  of  optimal  basis  sets  as  per  the 
Bagayoko,  Zhao,  and  Williams  (BZW)  procedure.17-20  As 
explained  elsewhere,  this  procedure  avoids  a  basis  set  and 
variational  effect  inherently  associated  with  variational  cal¬ 
culations  that  employ  a  basis  set  and  leads  to  the  calculated 
band  gaps  in  very  good  agreement  with  experiments.17-20 
With  the  above  method,  we  solved  the  Kohn-Sham14,15  equa¬ 
tion  self-consistently.  Self-consistency  was  followed  by  the 
calculations  of  the  total  energies.  Details  of  these  steps,  in¬ 
cluding  the  Kohn-Sham  equation  and  the  expression  for  the 
total  energy,  are  fully  described  in  the  literature.14-20 

In  the  LCAO  method,  we  expanded  the  electronic  eigen¬ 
function  'P*,'  of  the  many-atom  system  as  a  linear  combina¬ 
tion  of  atomic  wave  functions.17  These  input  functions  result 
from  ab  initio  calculations  for  atomic  species  that  are  present 
in  the  system.  For  the  calculations,  the  C(ls)  state  was  used 
as  the  core  state.  For  the  (10,  0)  and  the  (8,  4)  tubes,  the 
C(2s2p )  states  were  used  as  filled  and  partially  filled  valence 
states.  The  C(3s3p)  orbitals  were  the  unfilled  electron  states 
that  were  used  to  augment  the  basis  set  for  the  calculations. 
For  the  (10,  10)  tube,  the  C(2s2p)  states  were  used  as  filled 
and  partially  filled  valence  states  and  the  C(3s)  as  the  un¬ 
filled  state.  The  empty  C(3 p)  state  was  dropped,  in  the  case 
of  (10,  10),  due  to  convergence  difficulties. 

Uniaxial  strain  was  simulated  by  linearly  scaling  the  atom 
positions  along  the  tube  axis  in  the  carbon  nanotube.  To  find 
how  the  tube  radius  changes  with  uniaxial  strain,  after  the 
tube  axis  was  scaled,  the  radius  was  identified  from  the  total 
energy  minimization  procedure.  Namely,  for  a  given  nano- 
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TABLE  I.  'Unstrained  bond  lengths  (a0,  b0,  and  c0)  in  carbon 
nanotubes  (10,0),  (8,4),  and  (10,10).  L0  is  the  length  of  the  unit  cell; 
r0  is  the  tube  radius;  and  N  is  the  number  of  atoms  per  unit  cell. 


(10,0) 

(8,4) 

(10,10) 

a0  (A) 

1.416 

1.414 

1.420 

b0  (A) 

1.416 

1.419 

1.420 

cQ  (A) 

1.420 

1.420 

1.418 

r0  (A) 

3.915 

4.143 

6.781 

h0  (A) 

4.26 

11.27 

2.46 

N 

40 

112 

40 

tube  and  a  given  compression  or  tensile  strain,  the  total  en¬ 
ergy  was  obtained  as  a  function  of  the  tube  radius.  The  equi¬ 
librium  radius,  i.e.,  the  stable  one,  was  the  one  corresponding 
to  the  minimum  of  the  total  energy.  This  process  was  re¬ 
peated  for  each  strain  to  obtain  the  radius  of  a  tubule  as  a 
function  of  compression  or  tensile  strain.  Similarly,  the  cal¬ 
culated  bond  lengths,  band  gap,  and  Fermi  energies,  for  a 
given  nanotube  and  a  given  strain,  are  the  one  corresponding 
to  the  minimum  of  the  total  energies. 

The  Young’s  modulus  for  each  SWCNT  was  calculated  by 
fitting  the  total  energy/unit  cell  for  each  strain  ew  to  the 
equation 


E  = 


(1) 


where  Y  is  the  Young’s  modulus,  z  is  the  unit  cell  length,  and 
r  is  the  outer  radius  of  the  nanotube.  Equation  (1)  was  de¬ 
rived  using  the  assumption  that  the  nanotube  was  a  perfect 
cylinder.21  The  standard  radii  tabulated  for  the  nanotubes  are 
for  circumferences  through  the  centers  of  the  outer  atoms. 
The  outer  radius  includes  such  a  standard  radius  plus  the 
radius  of  the  carbon  atom  (0.71  A).  We  also  calculated  the 
in-plane  stiffness  C,  which  is  an  alternative  measure  of  the 
mechanical  characteristic  of  nanotubes  that  is  defined  as22 


C- 


1  &E 

~SdeV 


(2) 


where  S  is  the  surface  area  of  the  nanotube. 

Table  I  displays  the  unstrained  bond  lengths  and  other 
characteristics  of  the  tubules  under  study.  In  a  graphene 
sheet,  the  lengths  of  the  a,  b,  and  c  bonds  are  equal.  In 
unstrained  nanotubes,  which  are  rolled  graphene  sheets,  the 
three  bonds  are  generally  different.  For  the  SWCNT  (10,  0), 
the  c  bond  along  the  axis  is  greater  than  the  a  and  b  bonds 
along  the  circumference  ( a=b<c ).  For  the  SWCNT  (10, 
10),  the  c  bond  along  the  circumference  is  shorter  than  the  a 
and  b  bonds  along  the  axis  ( a=b>c ).  In  the  case  of  the 
SWCNT  (8,  4),  the  a  bond  closest  to  the  circumference  is 
less  than  the  c  bond  closest  to  the  axis,  i.e.,  a<b<c.  The 
bond  lengths  and  radii  as  functions  of  the  compression  or 
tensile  strain  are  displayed  in  Fig.  1  for  the  (10, 0)  and  Fig.  2 
for  the  (8,  4).  The  calculated  results  for  (10,  10)  are  very 
similar  to  those  of  (10,  0)  and  (8,  4)  and  are  not  presented  in 
a  separate  figure,  due  to  the  page  limitation  of  this  brief 
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FIG.  1 .  Bond  lengths  and  tube  radius  as  functions  of  compres¬ 
sion  and  tension  strain  e  for  (10,0)  nanotube.  Each  length  is  scaled 
to  its  unstrained  length. 

report.  As  per  Figs.  1  and  2,  the  calculated  and  normalized 
lengths  (i.e.,  a/aQ ,  h/h0,  or  c/cQ)  of  bonds  with  a  component 
along  the  axis  decrease  or  increase  linearly  with  compression 
or  tensile  strain,  respectively.  These  variations  are  quantified 
with  the  slopes  m  provided  in  Table  II. 

A  special  feature  of  the  calculated  properties  of  the  carbon 
nanotubes  was  that  the  radii  of  the  nanotubes  do  not  change 
under  a  substantially  large  uniaxial  strain,  from  -6%  to  +6%. 
The  Poisson’s  ratio,  which  is  the  ratio  of  the  transverse  con¬ 
tracting  strain  to  the  longitudinal  elongational  strain,  is 
nearly  zero  for  these  carbon  nanotubes.  This  behavior  is 
drastically  different  from  that  of  macroscopic  materials  and  it 
disagrees  with  results  of  tight-bonding  calculations  and  other 
mechanics  simulations.6,23  One  possible  explanation  is  that 
the  bonds  along  the  circumference,  partly  strengthened  by 
curvature  effects,  are  not  significantly  affected  by  uniaxial 
strains  orthogonal  to  them  (i.e.,  along  the  axis).  We  also  re¬ 
peated  the  studies  using  ab  initio  LDA  calculations  for 
SWCNT  (10,  0).  The  total  energy  minimization  of  the  LDA 
calculations  for  (10,  0)  found  the  bond  length  at  1.402  A  for 
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FIG,  2.  Bond  lengths  and  tube  radius  as  functions  of  tension 
strain  and  compression  strain  e  for  (8,4)  nanotube.  Each  length  is 
scaled  to  its  unstrained  length. 
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TABLE  II.  Coefficients  m  of  linear  variation  of  normalized  bond 
lengths  with  compression  (negative)  and  tensile  (positive)  strain 
(i.e.,  e  expressed  as  -0.04  and  +0.06  for  4%  and  6%  compression 
and  tensile  strain,  respectively).  Al/l0=me,  where  1  is  the  a,  b,  or 
c  bond  length.  The  Young’s  modulus  Y  is  in  TPa.  The  in-plane 
stiffness  C  is  in  J/m2. 


Nanotubes 

m  for  bond 

a 

m  for  bond 
b 

m  for  bond 

c 

Y 

(TPa) 

C 

(J/m2) 

(10,0) 

0.252 

0.252 

1.000 

1.47 

340 

(8,4) 

0.036 

0.572 

0.893 

1.10 

267 

(10,10) 

0.746 

0.746 

0a 

0.726 

272 

“Bond  c,  for  (10,10),  is  along  the  circumference  of  the  tube. 


a0  and  b0  and  1.406  A  for  c0,  that  is  about  1%  smaller  than 
the  results  of  the  ab  initio  GGA  calculations.  The  ab  initio 
LDA  calculations  for  SWCNT  (10,  0)  also  found  that  the 
radii  of  the  nanotube  did  not  change  under  a  substantially 
large  uniaxial  strain.  Among  the  previous  ab  initio  calcula¬ 
tions  of  the  elastic  properties  of  carbon  nanotubes,  Sanchez- 
Portal  et  al.  utilized  a  minimal  basis  set  of  one  s  and  three  p 
orbitals  per  carbon  atom  and  performed  LDA  calculations.24 
In  their  calculations,  they  used  pseudoatomic  orbitals.24  Van 
Lier  et  al.  utilized  ab  initio  Hartree-Fock  6-3 1G  method  and 
closed  nanotube  models  in  their  simulations.25  Sanchez- 
Portal  et  al.u  and  Van  Lier  et  al.2S  reported  relatively  small 
values  (from  0.14  to  0.19)  of  the  Poisson  ratio  for  their  cal¬ 
culated  carbon  nanotubes,  which  also  indicated  the  diameter 
rigidity  of  the  carbon  nanotubes.  These  previous  ab  initio 
calculations  utilized  different  computational  methods,  such 
as  a  minimal  basis  set  and  the  Hartree-Fock  method,  and 
their  results  are  slightly  different  from  ours.  As  demonstrated 
in  one  of  our  previous  publications,30  the  use  of  a  minimal 
basis  set  may  not  be  sufficient  to  obtain  a  highly  accurate 
solution  of  the  calculated  electronic  structure  of  the  carbon 
nanotubes,  which  could  explain  in  part  the  difference  be¬ 
tween  our  results  and  those  of  Sanchez-Portal  et  al  The  fun¬ 
damental  differences  between  the  Hartree-Fock  method  and 
density  functional  methods,  i.e.,  the  inclusion  of  the  elec¬ 
tronic  correlation  effects  in  the  latter,  partly  explain  the  dif¬ 
ference  between  our  results  and  those  of  Van  Lier  et  al.  In 
our  ab  initio  calculations,  we  utilized  extended  atomic  basis 
sets  and  performed  both  GGA  and  LDA  computations. 

The  calculated  Young’s  moduli  Y  (in  TPa)  for  the  tubes 
are  shown  in  Table  II.  We  recall  that  the  radius  included  in 
the  formula  for  the  Young’s  modulus  is  from  the  center  of  the 
tube  to  the  outer  circumference.  This  value  of  the  radius  led 
to  Young’s  moduli  close  to  the  experimentally  found  value  of 
approximately  1  TPa.9  The  (10,  0)  and  (10,  10)  tubes  have 
the  highest  and  lowest  Young’s  moduli,  respectively.  This 
trend  is  expected  due  to  the  bond  geometry  of  the  tubes 
according  to  chirality.  All  three  bonds  in  the  (10,  0)  tube 
have  a  significant  component  parallel  to  the  axis  along  which 
strain  is  applied.  In  the  case  of  (10,  0),  the  c  bond  is  entirely 
along  the  axis.  In  contrast,  one  of  the  bonds  of  the  (10,  10) 
tube  is  entirely  along  the  circumference;  it  is  not  expected  to 
oppose  any  resistance  neither  to  compression  nor  to  tensile 
strain.  For  this  reason,  nanotube  (10,  10)  is  intuitively  the 


£  (%) 

FIG.  3.  Band  gaps  of  indicated  carbon  nanotubes  as  functions  of 
compression  and  tensile  strain. 

easiest  to  stretch  or  to  compress,  as  confirmed  by  our  results 
in  Table  II. 

We  also  calculated  the  in-plane  stiffness  C  for  the  carbon 
nanotubes.  It  is  340,  267,  and  272  I/m2  for  (10,  0),  (8,  4), 
and  (10,  10),  respectively.  The  in-plane  stiffness  of  these 
nanotubes  exhibited  a  dependence  on  their  chirality,  when 
the  nanotubes  are  strained  along  the  tube  axis  (z-direction). 
Xiao  and  Liao  reported  an  average  in-plane  stiffness  C  of 
328  J/m2  for  graphene,  using  the  second-generation  Brenner 
potential  in  their  simulation  of  carbon  nanotubes.22 

The  calculated  electronic  properties  for  the  unstrained 
nanotubes  basically  reproduced  the  results  of  Zhao  et  al.20 
for  (10,  0)  and  (8,  4),  using  the  optimal  basis  sets  of  the 
BZW  method,  which  gave  converged  results  for  the  calcu¬ 
lated  electronic  structure  that  included  the  occupied  electron 
states  as  well  as  the  unoccupied  ones  near  the  Fermi  level. 
These  tubes  were  found  to  be  semiconductors.  These  authors 
provided  plots  of  the  band  structures  of  these  tubules.20  The 
electron  energy  bands  of  carbon  nanotube  (10,  10)  present  a 
semimetallic  property  as  reported  in  previous  publications. 
Figure  3  shows  the  nonlinear  variation  of  band  gaps  of  (10, 
0)  and  (8,  4)  with  strain.  For  both  nanotubes,  the  band  gap 
decreases  with  compression  strain  and  increases  with  tensile 
strain,  for  strain  values  smaller  than  or  equal  to  4%.  While 
the  band  gap  for  (8, 4)  reaches  a  minimum  for  a  compression 
strain  of  4%,  that  for  (10,  0)  exhibits  a  maximum  for  a  ten¬ 
sile  strain  of  4%.  Our  results  for  the  variation  of  the  band  gap 
with  strain  qualitatively  agree  with  the  findings  from  tight 
binding21  but  are  quantitatively  different,6  particularly  in  the 
magnitudes  of  the  band  gaps.  As  apparent  in  Fig.  3,  our 
results  indicate  that  nanotube  (10,  0)  become  metallic  at  6% 
compression  strain.  It  was  noticed  that  the  carbon  nanotubes 
may  collapse  for  large  strains,32,33  a  process  that  we  cannot 
simulate  at  present  using  the  ab  initio  quantum  calculations, 
because  the  required  computations  are  beyond  our  current 
computation  capability.  Therefore,  the  obtained  band  gaps 
for  large  strains  may  be  a  theoretical  simulation.  A  priori,  the 
general  behavior  of  these  band  gaps  with  strain,  particularly 
from  values  below  4%,  is  expected  on  the  basis  of  the  in¬ 
crease  or  decrease  of  the  overlap  between  atomic  sites,  for 
compression  or  tensile  strains,  respectively.  The  consequent 
broadening  or  flattening  of  the  bands,  respectively,  affects 
the  gaps. 

A  discussion  of  our  results  for  the  elastic  properties  is 
partly  limited  by  the  dearth  of  experimental  data  and  by  the 
rather  large  uncertainty  associated  with  currently  available 
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ones.  Indeed,  Treacy  et  al.26  and  Wong  and  coworkers27  re¬ 
ported  average  Young’s  moduli  of  1.8  and  1.28  TPa  with 
respective  uncertainties  of  1.4  and  0.59  TPa.  Given  these 
large  error  margins,  the  claimed  agreement  between  our  re¬ 
sults  and  experimental  ones  may  not  have  much  significance. 
A  comparison  with  previous  theoretical  studies  is  hampered 
by  the  differences  in  the  expressions  of  the  Young’s  modulus 
for  carbon  nanotubes  as  explained  by  Hernandez  and 
coworkers.10  These  authors  found  1.22  and  1.24  TPa  for  the 
Young’s  moduli  of  (10,  0)  and  (10,  10)  respectively.  The 
tight-binding  work  of  these  authors,  as  per  their  Fig.  3,  led  to 
about  the  same  modulus  for  (n , 0)  and  (n,n)  single-walled 
carbon  nanotubes  for  diameters  between  0.75  and  2  nm.  This 
finding  is  qualitatively  different  from  ours.  Lu9  also  found, 
using  the  empirical  force-constant  model,  values  of  0.975 
and  0.972  TPa  for  (10,  0)  and  (10,  10),  respectively.  Our 
results  are  clearly  different  for  the  three  nanotubes  under 
consideration.  Our  results  of  1.47  and  0.726  TPa  for  nano¬ 
tubes  (10,  0)  and  (10,  10),  respectively,  qualitatively  agree 
with  the  trend  in  the  recent  tight-binding  work  of  Zhang  et 
al.2S  who  found  the  elastic  limit  of  (n,0)  tubes  to  be  about 
twice  that  of  ( n,n )  tubes  of  comparable  radius.  Another  in¬ 
dication  of  the  possible  correctness  of  our  findings  may  re¬ 
side  in  the  case  of  graphite.  Indeed,  the  Young’s  modulus  of 
graphite  along  the  c  axis,  0.0365  TPa,  is  very  different  from 
the  one  in  the  basal  plane,  1.02  TPa.29 

Due  to  our  utilization  of  the  BZW  method,  our  calculated 
band  gaps  are  expectedly  higher  than  other  theoretical  find¬ 
ings  known  to  us.  While  Reich  et  al .30  reported  a  gap  of 
0.8  eV  for  nanotube  (8,  4),  we  found  0.96  eV.  More  impor¬ 
tantly,  our  results  clearly  show  that  the  variation  of  the  band 
gap  with  strain  is  far  from  being  linear.  For  nanotubes  (10,  0) 
and  (8,  4),  the  average  increase  of  the  band  gap  per  tensile 
strain  under  2%  is  0.125  eV.  This  result  is  larger  than  a  pre¬ 
viously  estimated  maximum  variation31  of  0.100  eV  per  per¬ 
cent  strain,  as  is  reportedly  the  case  for  bulk  semiconductor. 

In  conclusion,  this  work  reported  the  calculated  structural, 
elastic,  and  electronic  properties  of  selected,  single-wall  car¬ 
bon  nanotubes  (10,  0),  (8,  4),  and  (10,  10)  under  uniaxial 


strain.  The  results  of  our  ab  initio,  self-consistent,  GGA- 
BZW  calculations  agreed  qualitatively  with  some  findings  of 
the  tight  binding  approach  and  quantitatively  disagreed  with 
them.  We  found  that  the  radii  of  these  nanotubes  do  not 
change  under  uniaxial  strain,  up  to  6%.  We  quantified  the 
linear  variation  of  bond  lengths  with  compression  and  tensile 
strain  and  showed  the  change  of  the  band  gap  with  strain,  for 
semiconducting  nanotubes,  to  be  essentially  nonlinear.  Mea¬ 
surements  of  elastic  limits  for  (n,0)  and  (n , n)  nanotubes  and 
of  the  Young’s  moduli  and  the  in-plane  stiffness  for  graphite 
qualitatively  agree  with  the  trends  in  our  calculated  results 
for  the  selected  single-walled  carbon  nanotubes. 

Recently,  we  received  two  preprints  regarding  the  experi¬ 
mental  work  by  S.  B.  Cronin  et  al.  at  Harvard  University  and 
the  Massachusetts  Institute  of  Technology  on  measurements 
of  uniaxial  strain  in  single-wall  carbon  nanotubes,  utilizing 
resonance  Raman  spectra  of  atomic-force  microscope  modi¬ 
fied  SWCNTs.  This  work  entailed  bending  the  nanotubes 
while  holding  the  ends  fixed.  Two  important  comments 
should  be  made  about  this  finding.  The  first  one  consists  of 
the  fact  that  our  calculations  did  not  include  a  determination 
of  whether  or  not  the  nanotubes  actually  conserve  their  sym¬ 
metrical  geometry  under  strain  as  high  as  6%.  Cronin  et  al.32 
reported  cases  where  carbon  nanotubes  broke  under  strains 
greater  than  1.65%.  The  second  comment  stems  from  an  in¬ 
dication,  from  the  work  of  Cronin  et  al. , 32,33  that  the  radii  of 
semiconducting  SWCNTs  do  not  change  under  strains  be¬ 
tween  0.06  and  1.65%.  The  constancy  of  the  radial  breathing 
mode  (RBM)  frequency  (wRBM)  lead  to  this  conclusion,32,33 
given  that  the  tube  diameter  d,=248/  o>RBM. 

ACKNOWLEDGMENTS 

This  work  was  funded  in  part  by  the  National  Aeronautics 
and  Space  Administration  (NASA  Award  No.  NCC  2-1344), 
the  Department  of  the  Navy,  Office  of  Naval  Research 
(ONR,  Awards  No.  N00014-05-1-0009  and  No.  N00014-04- 
0587),  and  the  National  Science  Foundation  (NSF  Award  No. 
HRD-0000272). 


*Corresponding  author.  Email  address:  zhao@grant.phys.subr.edu 
'M.  Wilson,  K.  Kannangara,  G.  Smith,  M.  Simmons,  and  B.  Ra- 
guse,  Nanotechnology:  Basic  Science  and  Emerging  Technolo¬ 
gies  (Chapman  and  Hall,  Boca  Raton,  2002). 

2C.  T.  White  and  J.  Mintmire,  Nature  (London)  394,  29  (1998). 
3J.  W.  Mintmire  and  C.  T.  White,  Phys.  Rev.  Lett.  81,  2506 
(1998). 

4L.  Yang,  M.  P.  Anantram,  J.  Han,  and  J.  P.  Lu,  Phys.  Rev.  B  60, 
13  874  (1999). 

5L.  Yang  and  J.  Han,  Phys.  Rev.  Lett.  85,  154  (2000). 

6L.  Yang,  J.  Han,  M.  P.  Anantram,  and  R.  Jaffe,  Comput.  Model. 
Eng.  Sci.  3,  675  (2002). 

7M.  Ouyang,  J.  L.  Huang,  C.  L.  Cheung,  and  C.  M.  Lieber, 
Science  292,  702  (2001). 

8B.  I.  Yakobson,  C.  J.  Brabec,  and  J.  Bemholc,  Phys.  Rev.  Lett. 
76,  2511  (1996). 


9J.  P.  Lu,  Phys.  Rev.  Lett.  79,  1297  (1997). 

I0E.  HernSndez,  C.  Goze,  P.  Bernier,  and  A.  Rubio,  Phys,  Rev.  Lett. 
80,  4502  (1998). 

11 D.  Sanchez-Portal,  E.  Artaeho,  J.  M.  Soler,  A.  Rubio,  and  P. 

Ordejdn,  Phys.  Rev.  B  59,  12  678  (1999). 

12V.  N.  Popov,  V.  E.  Van  Doren,  and  M.  Balkanski,  Phys.  Rev.  B 
61,  3078  (2000). 

13J.  P.  Perdew  and  W.  Yue,  Phys.  Rev.  B  33,  8800  (1986);  J.  P. 
Perdew,  ibid.  33,  8822  (1986);  J.  P.  Perdew  and  A.  Zunger,  ibid. 
23,  5048  (1981). 

14P.  Hohenberg  and  W.  Kohn,  Phys.  Rev.  136,  B864  (1964). 

15W.  Kohn  and  L.  J.  Sham,  Phys.  Rev.  140,  A1133  (1965). 

16 J.  Callaway  and  N.  H.  March,  in  Solid  State  Physics,  Vol.  38, 
edited  by  H.  Ehrenreich,  D.  Turnbull,  and  F.  Seitz  (Academic 
Press,  New  York,  1984),  p.  135,  and  references  therein. 
i7G.  L.  Zhao,  D.  Bagayoko,  and  T.  D.  Williams,  Phys.  Rev.  B  60, 


205410-4 


STRUCTURAL,  ELASTIC,  AND  ELECTRONIC... 

1563  (1999). 

18D.  Bagayoko,  G.  L.  Zhao,  J.  D.  Fan,  and  J.  T.  Wang,  J.  Phys.: 

Condens.  Matter  10(25),  5645  (1998). 

19D.  Bagayoko  and  G.  L.  Zhao,  Int.  J.  Mod.  Phys.  B  13(29,  30  & 
31),  3767  (1999). 

20G.  L.  Zhao,  D.  Bagayoko,  and  L.  Yang,  Phys.  Rev.  B  69,  245416 
(2004). 

21 M.  Dresselhaus,  G.  Dresselhaus,  and  P.  Ecklund,  Science  of 
Fullerenes  and  Carbon  Nanotubes  (Academic  Press,  San  Diego, 
1996);  also  see  C.  Kittel,  Introduction  to  Solid  State  Physics ,  7th 
ed.  (Wiley,  Brisbane,  1996). 

22T.  Xiao  and  K.  Liao,  Phys.  Rev.  B  66,  153407  (2002). 

23  A.  Sears  and  R.  C.  Batra,  Phys.  Rev.  B  69,  235406  (2004). 

24D.  Sanchez-Portal,  E.  Artacho,  J.  M.  Soler,  A.  Rubio,  and  P. 
Ordejon,  Phys.  Rev.  B  59,  12  678  (1999). 

25  G.  Van  Lier,  C.  Van  Alsenoy,  V.  Van  Doren,  and  P.  Geerlings, 

Chem.  Phys.  Lett.  326,  181  (2000). 

26  M.  M.  J.  Treacy,  T.  W.  Ebbesen,  and  J.  M.  Gibson,  Nature 


PHYSICAL  REVIEW  B  71,  205410  (2005) 
(London)  381,  678  (1996). 

27E.  W.  Wong,  P.  E.  Sheehan,  and  C.  M.  Lieber,  Science  277,  1971 

(1997). 

28Peihong  Zhang,  Paul  E.  Lammert,  and  Vincent  H.  Crespi,  Phys. 
Rev.  Lett.  81,  5346  (1998). 

29 O.  L.  Blakslee  etal.,  J.  Appl.  Phys.  41,  3373  (1970);  E.  J.  Seldin 
and  C.  W.  Nezbeda,  ibid.  41,  3389  (1970). 

30S.  Reich,  C.  Thomsen,  and  P.  Ordejon,  Phys.  Rev.  B  65,  155411 

(2002). 

31 E.  D.  Minot,  Yuval  Yaish,  Vera  Sazonova,  Ji-Yong  Park,  Markus 
Brink,  and  Paul  L.  McEuen,  Phys.  Rev.  Lett.  90(15),  156401 
(2003). 

32S.  B.  Cronin,  A.  K.  Swan,  M.  S.  Unlu,  B.  B.  Goldberg,  M.  S. 
Dresselhaus,  and  M.  Tinkham,  Phys.  Rev.  Lett.  93,  167401 
(2004). 

33 S.  B.  Cronin,  A.  K.  Swan,  M.  S.  Unlu,  B.  B.  Goldberg,  M.  S. 
Dresselhaus,  and  M.  Tinkham  (unpublished). 


205410-5 


Ab-Initio  Simulations  of  the  Growth  and  Structural  Properties  of  Short 

Carbon  Nanobells 

G.  L.  Zhao1,a,  D.  Bagayoko1,  and  E.  G.  Wang2 

’Department  of  Physics,  Southern  University  and  A  &  M  College 
Baton  Rouge,  Louisiana  70813,  USA 
a  E-mail  address:  zhao@grant.phys.subr.edu 
2State  Key  Laboratory  for  Surface  Physics,  Institute  of  Physics, 

Chinese  Academy  of  Sciences,  Beijing,  100080,  P.  R.  China 

Keywords:  carbon  nanobells,  carbon  nanotubes,  growth  and  structural  properties  of  nanobells. 

Abstract:  We  performed  ab-initio  density  functional  simulations  to  study  the  structural  and 
growth  properties  of  short  carbon  nanobells.  We  used  a  real  space  approach  and  the  linear 
combination  of  atomic  orbitals  (LCAO)  formalism.  In  the  nitrogen-doped  carbon  nanobells,  the 
nitrogen  atoms  that  are  attracted  to  the  open-edge  sites  of  the  carbon  nanobells  play  an  important 
role  in  the  growth  of  the  short  carbon  nanostructures.  We  present  the  calculated  electronic 
structure  of  the  short  nanobells.  The  calculated  local  density  of  states  of  the  nanobells  revealed 
field  emission  characteristics  that  agree  with  experimental  observations. 


I.  Introduction 

The  carbon  nanostructures  present  interesting 
properties  with  a  great  potential  of 
applications.  Particularly,  short  carbon 
nanotubes  or  carbon  nanobells  present 
excellent  field-emission  properties  that  have 
attracted  a  great  interest  for  applications.1'8 
Experiments  demonstrated  that  the  nitrogen- 
assisted  synthesis  can  grow  carbon  nanofibers 
on  a  large  scale.9'12  Such  carbon  nanofibers 
exhibit  a  "bamboo-like"  structure. 
Distinctively,  a  great  part  of  the  bamboo-like 
nanofibers  consists  of  short  carbon  nanobells. 
Individual  nanobells  are  self-contained  and 
stack  one  on  top  of  the  other  to  create  a  long 
nanofiber.  The  carbon  nanobells  may  be 
viewed  as  short  nanotubes  such  that  their 
lengths  are  of  the  same  order  as  their 
diameters.  The  electronic  structure  of  the 
short  carbon  nanobells  is  substantially 
different  from  that  of  the  pure  and  long 
carbon  nanotubes  because  of  the  quantum 
effects  inherent  to  their  sizes.  The  nanobells 
exhibit  novel  electron  field  emission 


properties  with  a  turn-on  field  of  electron 
emission  as  low  as  0.8  V  /  jum  ' 9,13  We 
performed  ab-initio  density  functional 
calculations,  aiming  to  understand  the 
growth,  structural,  and  electronic  properties 
of  short  carbon  nanobells. 

II.  Method 

Although  the  experimental  synthesis  of 
carbon  nanobells  produced  samples  of 
various  sizes,  measurements  revealed  that  the 
atomic  structure  of  the  bamboo-like 
morphology  in  thick  nanobells  is  similar  to 
that  of  the  thinner  ones.9  Experiments  also 
revealed  that  adjacent  nanobells  do  not  have  a 
firm  contact,  but  instead  the  closed  end  of  one 
nanobell  is  weakly  inserted  into  the  open  end 
of  another.  Single  nanobells  can  be  easily 
separated  from  other  part  of  the  nanofiber.10 
Hence,  our  ab-initio  calculations  will  focus 
on  single  nanobells.  The  weak  interaction 
between  adjacent  carbon  nanobells  may  be 
included  as  a  perturbation  in  further  studies  in 
modeling  the  bamboo-like  morphology  of  the 
nanofibers.  We  constructed  a  prototypical 
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carbon  nanobell  as  a  model  to  simulate  the 
general  structure  of  short  carbon  nanobells. 
Fig.  1  presents  the  calculated  images  of  the 
electron  density  distribution  of  a  prototypical 
carbon  nanobell  in  real  space.  The 
prototypical  carbon  nanobell  includes  150 
atoms  and  consists  of  half  of  the  C240 

fullerene  with  30  extra  atoms  on  the  bell.  Its 

0 

length  and  diameter  are  8.5  and  7.4  A, 
respectively.  Although  the  size  of  this 
prototypical  carbon  nanobell  may  differ  from 
that  of  fabricated  samples,  their  fundamental 
features  are  similar.  These  features  include 
the  structural  properties  of  the  open-edge,  the 
wall,  and  the  closed  cap  of  the  nanobells. 


Fig.  1.  The  calculated  images  of  the  electron 
density  distribution  in  real  space  for  a 
prototypical  carbon  nanobell  that  is  viewed  from 
different  directions 


Our  ab-initio  calculations  are  based  on 
the  density  functional  theory  of  Hohenberg- 
Kohn  and  Kohn-Sham.14, 15,16,17,18  In  the 
linear  combination  of  atomic  orbitals  (LCAO) 
method,  we  solved  the  Kohn-Sham  equations 
self-consistently  by  employing  the  Rayleigh- 
Ritz  variational  process.19,20  In  the  self- 
consistent  calculations  for  the  electronic 
structure  of  short  CNX  nanobells,  we 
employed  an  extended  basis  set  that  includes 
atomic  orbitals  of  C(ls2s3s  2p3p)  and 
N(ls2s3s  2p3p).  Here  C(3s  3p)  and  N(3s  3p) 
are  extra  orbitals  that  are  used  to  augment  the 
basis  set  to  account  for  possible  charge 
diffusion  and  polarization  in  the  short 
nanobells.  The  real  space  approach  of  the 
LCAO  method  enables  us  to  complete  the 
required  computations  using  our 
SiliconGraphics  0rigin2000  that  is  equipped 
with  2  GB  RAM  (memory). 


III.  Results 

(A).  Relaxed  Atomic  Structure  near  the 
Open-Edge  of  the  Short  Carbon  Nanobell. 

The  open-edges  of  carbon  nanobells  play  an 
important  role  in  their  electronic  properties. 
The  structure  near  the  open  edge  of  the 
carbon  nanobell  is  different  from  that  in  the 
inner  wall  and  the  closed  cap.  We  performed 
total  energy  minimization  to  identify  the 
atomic  structure  near  the  open-edge  of  the 
prototypical  carbon  nanobell.  The  calculated 
C-C  bond  length  on  the  first  atomic  ring  at 
the  open-edge  of  the  carbon  nanobell  is 

o 

1.37 A,  which  is  much  shorter  than  the  C-C 

bond  length  of  1.415  A  in  the  inner  atomic 
rings  of  the  bell.  This  reduction  of  the  C-C 
bond  length  for  the  carbon  atoms  at  the  open- 
edge  sites  is  a  result  of  the  reduction  in 
coordination  number.  The  calculated  bond 
length  of  the  carbon  atoms  on  the  first  atomic 
ring  to  the  carbon  atoms  on  the  second  ring  is 

o 

1.40  A ,  which  is  not  much  different  from  that 
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in  the  inner  rings,  away  from  the  open-edge 
of  the  nanobell.  This  result  indicates  that  the 
effect  of  further  relaxation  of  the  third  atomic 
ring  and  inner  rings  would  be  insignificant. 

(B).  A  Possible  Growth  Mechanism  of 
Carbon  Nanobells.  The  growth  mechanism 
of  the  carbon  nanobells  remains  a  difficult 
problem.  There  are  probably  two  major 
growth  models.1  The  first  one  assumes  that 
the  carbon  atoms  are  added  at  the  open-ends 
of  the  nanobells.  1,11  The  second  one  involves 
the  C2  absorption  process  that  is  assisted  by 
the  pentagonal  defects  on  the  closed  caps.1 
Although  both  models  are  very  interesting, 
we  have  technical  difficulties  using  density 
functional  computations  to  simulate  the 
second  growth  model.  We  successfully 
performed  local  density  functional 
computations  to  test  the  first  growth 
mechanism.  Because  of  the  active  dangling 
bonds  of  carbon  atoms  at  the  open-edge  of 
nanobells,  there  is  a  high  possibility  that 
carbon  atoms  can  be  attracted  to  these  sites 
for  growth  to  occur.  We  calculated  the  total 
energies  in  the  following  two  cases.  In  the 
first  case,  a  ring  of  20  carbon  atoms  attaches 
to  the  open-edge  of  the  carbon  nanobell  at  a 

o 

bond-length  of  1.415  A.  In  the  second  case, 
these  20  carbon  atoms  are  free.  The 
calculated  total  energy  in  case  I  is 
substantially  lower  than  the  corresponding 
value  in  case  II.  The  total  energy  difference, 
which  may  also  be  defined  as  the  cohesive 
energy  of  the  carbon  atoms  in  the  first  growth 
model,  is  4.7  eV/atom.  In  this  growth  model, 
see  Fig.  1,  C2  dimers  are  absorbed  at  the 
active  dangling  bond  sites  at  the  open  end  of 
the  carbon  nanobell.  A  C2  dimer  that  deposits 
on  the  open-edge  of  the  bell  forms  one 
covalent  bond  on  the  same  atomic  ring.  Each 
of  the  C  atoms  also  forms  one  covalent  bond 
with  the  C  atom  of  the  next  atomic  ring  and 
remains  an  active  dangling  bond  towards  the 
open-space.  This  grow  model  is  consistent 
with  previous  analysis  of  experimental 


results. n’13,  21  Although  we  cannot  rule  out 
the  second  growth  model  (growth  on  the 
closed  cap  of  carbon  nanobells),  our 
calculations  presented  a  clear  evidence  that 
the  first  growth  mechanism  is  highly  possible. 

(C).  Nitrogen-Doped  Carbon  Nanobell. 

Experimental  results  indicated  that  the  growth 
of  the  short  carbon  nanobells  is  highly 
dependent  on  the  nitrogen  concentration  in 
the  gas  mixture  during  the  synthesis.9,10 
Without  nitrogen  atoms  in  the  growth  gas 
mixture  under  the  same  conditions,  long  and 
pure  carbon  nanotubes,  without  the  “bamboo¬ 
like”  morphology,  are  produced.  However,  it 
was  not  clear  why  the  nitrogen  atoms  could 
turn  the  growth  of  the  would-be  long  carbon 
nanotubes  into  that  of  short  carbon  nanobells. 
We  performed  total  energy  calculations  to 
study  the  nitrogen-doped  carbon  nanobells 
(CNX  nanobells).  Since  nitrogen  atoms  can 
also  form  the  planar  sp  -hybrid,  we  studied 
the  substitutional  doping  of  nitrogen  atoms  in 
the  CNX  nanobells.  The  substitutional  doping 
of  nitrogen  atoms  in  CNX  nanobells  was  also 
proposed  from  the  analysis  of  experimental 
results.13 

We  first  replaced  carbon  atoms  with 
nitrogen  atoms  at  the  open-edge  of  the 
prototypical  nanobell.  We  compared  the  total 
energies  in  the  two  cases.  In  the  first  case, 
computations  are  carried  out  for  the  pristine 
carbon  nanobell  in  the  presence  of  ten  free 
atoms  of  nitrogen.  In  the  second  case,  ten 
nitrogen  atoms  replace  ten  carbon  atoms  on 
the  open-edge  and  the  substitutionally  N- 
doped  nanobell  is  in  the  presence  of  ten  free 
atoms  of  carbon.  The  calculated  total  energy 
in  the  second  case  is  lower  than  that  in  the 
first  case  by  an  amount  of  0.5  eV  per  atom, 
without  relaxing  the  structure.  We  then 
performed  total  energy  calculations  and 
identified  the  atomic  positions  of  the  relaxed 
structure  of  the  nitrogen-doped  nanobell  from 
the  total  energy  minimization.  The  calculated 
C-N  bond  length  in  the  first  atomic  ring  at  the 


3 


open-edge  of  the  nanobell  is  1.376  A.  The 
bond-length,  from  the  nitrogen  atoms  of  the 
first  atomic  ring  to  the  carbon  atoms  of  the 

o 

second  atomic  ring,  is  1.377  A.  They  are 
close  to  the  length  of  the  partially  double  C-N 

o 

bond  (1.352 A)  in  heterocyclic  systems. 
The  C-C  bond  length  from  the  first  atomic 

o 

ring  to  the  second  ring  is  1.397  A .  The  total 
energy  of  the  nitrogen-doped  nanobell,  in  the 
relaxed  structure,  is  lower  than  that  of  the 
undoped  one  by  0.72  eV  per  atom. 

We  also  studied  the  substitutional 
doping  with  nitrogen  atoms  (for  carbon 
atoms)  in  the  wall  (away  from  the  open-edge 
and  the  closed  cap)  of  the  prototypical  carbon 
nanobell.  We  used  the  same  computation 
procedure  as  discussed  above.  The  calculated 
total  energy  of  the  N-doped  nanobell,  at  the 
wall  position,  is  much  higher  than  that  of  the 
undoped  bell  by  an  amount  of  5  eV  per  atom. 
This  result  indicates  that  the  substitutional 
doping  with  nitrogen  in  the  wall  of  carbon 
nanobells  is  not  energetically  favored. 

These  calculated  results  for  the 
nitrogen-doped  nanobells  in  the  prototypical 
model  thus  indicate  that  the  dopant  nitrogen 
atoms  prefer  to  stay  at  the  open  edge  of  the 
nanobell  as  opposed  to  being  in  the  wall  (i.e. 
lower  rings)  of  the  nanobell.  Once  there  are 
enough  nitrogen  atoms  on  the  open  edge  of 
the  nanobell,  carbon  atoms  cannot  attach  to 
the  open-edge  of  the  CNX  nanobell.  Such  an 
attachment  will  result  in  nitrogen  atoms  being 
in  inner  rings  other  than  that  at  the  open-edge, 
resulting  in  a  configuration  that  is  not  favored 
energetically.  Consequently,  nitrogen  atoms 
act  as  the  stopper  of  the  growth  of  the  carbon 
nanobell.  This  growth  mechanism  for  the 
nitrogen-doped  carbon  nanobells  and  the 
effect  of  the  nitrogen  atoms  in  the  growth 
process  agree  with  the  observations  from 
experiments.9)10,1 1,13 
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Fig.  2.  The  calculated  local  density  of  states  of 
electrons  of  a  nitrogen-doped  carbon  nanobell. 
The  1st  ring-N  and  -C  refer  to  the  nitrogen  and 
carbon  atoms  on  the  first  atomic  ring  at  the  open- 
edge  of  the  nanobell.  The  2nd  and  3rd  rings  refer 
to  the  second  and  third  atomic  rings  near  the 
open-edge.  The  inner  rings  refer  to  the  atomic 
sites  away  from  the  open-edge,  including  the 
closed  cap.  The  Fermi  level  is  at  0.0  eV. 


We  present  the  calculated  local 
density  of  states  (LDOS)  of  the  CNX  nanobell 
in  Fig.  2.  The  dopant  nitrogen  atoms  stay  at 
the  open-edge  sites  in  the  prototypical 
nanobell.  The  LDOS  from  the  nitrogen  and 
carbon  atoms  on  the  first  atomic  ring  at  the 
open-edge  of  the  nanobell  are  shown  in  the 
first  and  the  second  panels  of  Fig.  2, 
respectively.  There  is  no  N  atom  on  other 
atomic  rings  away  from  the  open-edge  of  the 
nanobell.  Fig.  2  clearly  shows  that  the  N  and 
C  atoms  on  the  first  atomic  ring  at  the  open- 
edge  of  the  nanobell  have  the  dominant 
contribution  to  the  electronic  states  around 
the  Fermi  level.  The  LDOS  at  the  Fermi 
level,  from  the  N  and  C  atoms  on  the  first 
atomic  ring  at  the  open-edge,  are  0.56  and 
0.71  states  per  eV  per  atom,  respectively. 
The  contribution,  from  the  second  atomic  ring 
near  the  open-edge,  to  the  LDOS  at  the  Fermi 
level  is  0.23  states  per  eV  per  atom.  The  third 
atomic  ring  and  the  inner  rings,  that  include 
the  closed  cap,  have  a  much  smaller 
contribution  to  LDOS  at  the  Fermi  level,  at 
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about  0.08  states  per  eV  per  atom.  Thus, 
electrons  will  be  emitted  from  the  atomic 
sites  near  the  open-edge  of  the  nanobell  in 
field-emission  experiments.  The  LDOS  of 
Fig.  2  exhibit  a  nearly  metallic  behavior.  The 
real  nanobell  samples  in  experimental  studies 
mostly  involve  multilayers  of  graphite  sheets 
and  the  nanobells  stacked  one  on  top  of  the 
other.  The  weak  interactions  between  graphite 
layers  in  these  nanobells  and  that  between  the 
bells  may  further  broaden  the  electronic 
energy  levels.  Consequently,  a  metallic 
behavior  may  be  observable  in  the 
measurements  of  the  electronic  structure  of 
N-doped  nanobells. 

IV.  Conclusion 

In  conclusion,  we  performed  ab-initio 
density  functional  calculations  to  simulate  the 
growth,  structural,  and  electronic  properties 
of  prototypical  carbon  nanobells.  In  the 
nitrogen-doped  carbon  nanobells,  nitrogen 
atoms  that  are  attracted  to  the  open-edge  of 
the  nanobells  may  play  a  role  of  stopper  of 
the  growth  of  the  nanostructures.  The 
calculated  local  densities  of  states  of  the  CNX 
nanobells  indicate  that  electrons  are  most 
likely  emitted  from  the  atomic  sites  near  the 
open-edge  of  the  carbon  nanobell  in  field- 
emission  experiments.  This  result  agrees  with 
experimental  observations. 
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Abstract.  We  present  a  universal  relation  between  the  densities  of  states  near  Van  Hove 
singularities  and  the  effective  electron  masses  in  1 -dimensional  (1-D)  semiconductors. 
The  relation  can  be  utilized  as  a  new  method  to  determine  the  effective  masses  of  charge 
carriers  of  1-D  semiconductors  in  theoretical  calculations  and  in  experimental 
measurements.  The  calculated  results,  utilizing  the  relation  for  1-D  single-walled  carbon 
nanotubes  (SWCNTs),  agree  well  with  that  of  the  conventional  calculations  when  both 
approaches  utilize  the  outputs  of  ab-initio  density  functional  computations. 


I.  Introduction. 

One-dimensional  (1-D)  materials,  such 
as  carbon  nanotubes  (CNTs),  have 
emerged  as  attractive  materials  for 
applications  in  molecular-scale 
electronics.1,11  The  electronic  properties 
of  1-D  CNTs  have  been  extensively 
studied  in  recent  years  using  various 
experimental  tools  and  theoretical 
calculations.11,111,lv  Particularly,  scanning 
tunneling  microscopy  (STM)  is  a 
favorable  tool  for  measuring  atomic 
structures,  density  of  states  (DOS)  of 
electrons,  Van  Hove  singularities,  and 
energy  band  gaps.  However,  there  are 
some  other  physical  quantities  that 
require  further  studies.  One  of  them  is 
the  effective  mass  of  charge  carriers 
(electrons  or  holes),  an  important 
physical  quantity  that  characterizes  the 
electronic  transport  properties  of 
materials.  Here,  we  present  a  universal 
relation  between  the  densities  of  states 
of  electrons  near  Van  Hove  singularities 
and  effective  electron  masses  in  1-D 


materials.  This  relation  can  be  used  as  a 
new  method  to  determine  the  effective 
masses  of  charge  carriers  in  any  1-D 
semiconductors. 


II.  Densities  of  States  and  Effective 
Masses  in  1-D  Materials. 

The  electron  density  of  states  (DOS)  per 
unit  cell,  G(b),  of  a  1-D  material  can  be 
calculated  by  the  following  expression/ 


G(e)  = 


—  X 

O  TT 

Z n  (%=«) 


1 

fctk 

dk 


(1) 


where  ejk  is  the  energy  of  the  i-th 
electron  band  at  point  k  in  the 
irreducible  Brillouin  zone;  and  L  is  the 
length  of  the  unit-cell. 

At  an  extremal  point  of  an 
electron  energy  band,  the  first  derivative 
of  the  band  energy  with  respect  to  the 
wave  vector  k  is  zero,  that  is  deik  /dk  = 

0.  A  Van  Hove  singularity  occurs  at 
this  point.  We  can  then  approximate  the 
band  dispersion  as  a  parabolic 


l 


expression  in  the  neighborhood  of  the 
extremal  point  that  particularly  includes 
the  top  or  the  bottom  of  an  energy  band, 
£ik  =£i±ai(k-ki)2  (2) 


where  k,  and  s,  are  the  location  and 
the  band  energy  at  the  extremal  point, 
respectively.  Here  a,  >  0.  When  the 
second  derivative  of  the  band  energy 
with  respect  to  the  wave  vector  k  is  also 
zero,  the  effective  mass  of  the  charge 
carrier  near  the  Van  Hove  singularity  is 
not  defined.  Therefore,  we  do  not 
discuss  this  case  further  in  this  article. 

We  can  then  calculate  the 
electron  density  of  states  per  unit  cell 
near  the  Van  Hove  singularity  as 
follows 

D(e)  =  G(e)/L  = 


(3) 


where  D(e)  is  a  reduced  density  of  states, 
which  does  not  depend  on  the  length  of 
the  unit  cell. 

We  can  calculate  the  effective 
mass  of  charge  carriers  using  the 
electron  energy  band  dispersion, 


h2 
i2  £ 

1  cik 

dk2 


h2 

2a, 


(4) 


where  m*  is  the  effective  mass  of  charge 


carriers  (electrons  or  holes)  near  the 
extremal  point  of  the  i-th  band. 


band  as  D](e)  = 


4 n 


and 


m!  = 


2  ax 


If  a  second  single  band  of  the  same 
material  or  another  1-D  material  near  the 
Van  Hove  singularity  is, 

£ Ik  ~  £"l  —  ^2  —  ^2  )  (6) 

We  then  have  the  density  of  states  and 
the  effective  mass  of  this  second  single 

band  as  D2(s)  =  * 


m2  = 


r 

2a , 


and 


Using  the  above  equations,  we 
identify  the  following  relation, 


m, 


w, 


a  i 


£>i  (g  +  g0 ) 

DJ£  +  e2) 


(7) 


The  above  relation  does  not 
depend  on  the  energy.  Once  the  effective 
mass  for  a  single  electron  band  of  a  1-D 
material  is  known,  we  can  then  take  it  as 
a  reference.  The  effective  masses  of 
charge  carriers  (electrons  or  holes)  for 
any  other  single  bands  of  other  1-D 
materials  can  be  calculated  using  the 
above  relation  and  corresponding 
densities  of  states  near  Van  Hove 
singularities. 


III.  A  Relation  between  DOS  and  m* 
in  1-D  Materials. 

(A).  Case  of  single  bands.  We  first 
consider  the  case  of  a  single  band.  The 
dispersion  of  a  single  electron  band  of  a 
1-D  material  near  its  Van  Hove 
singularity  can  be  written  as 

£lk  =£,±al(k-k1)2  (5) 

We  obtain  the  density  of  states  and  the 
effective  electron  mass  of  this  single 


(B).  Case  of  bands  that  are  degenerate 
at  an  extremal  point,  but  are  not 
degenerate  elsewhere.  If  two  bands  are 
degenerate  at  an  extremal  k-point,  but 
are  not  degenerate  elsewhere,  we  may 
assume  the  dispersions  of  the  two  as, 
£f=£2±fr(0(*-*2)2  (8) 

where  i=l  or  2  for  the  two  bands, 
respectively.  We  may  take  a  single  band 
of  a  1-D  material  as  a  reference,  for 
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which  the  band  dispersion  follows  Eq. 
(5).  The  density  of  states  and  the 
effective  mass  of  this  single  band  has 
been  discussed  in  the  previous  section. 
We  then  identify  the  relation  as, 

•P2(g  +  gi)  /m 

2D1(e  +  £0)  _ 


HI, 


where  m*  and  D2(e)  are  the  average 
effective  mass  and  the  density  of  states 
of  the  doubly  degenerate  bands, 
respectively. 


(10) 


where  mf  = - and  mi2)  = - . 

2  2b,  2  2  b2 

The  above  discussions  can  also 
be  extended  to  three  or  more  energy 
bands  that  are  degenerate  at  an  extremal 
point,  but  are  not  degenerate  elsewhere. 


IV.  Applications. 

Although  there  is  a  wealth  of 
experimental  results  regarding  the 
density  of  states  of  electrons  in  single- 
walled  carbon  nanotubes  (SWCNTs), U,U1 
the  study  of  the  effective  masses  of  the 
charge  carriers  in  these  materials  is  far 
from  complete.  As  a  test  case  of  the  new 
method  to  study  the  effective  masses  of 
the  charge  carriers  in  1-D  materials,  we 
applied  Eq.(7)  or  (9)  to  some  SWCNTs. 
We  list  some  of  the  results  in  Table  1, 
where  mo  is  the  free  electron  mass.  Ef’ 
in  Table  1  refers  to  the  first  valence  band 
near  the  Fermi  level  (Ef),  whereas  E’f 

refers  to  the  first  conduction  bands  from 
Ef.  Em  is  the  energy  (in  eV)  of  the  top 
of  a  valence  band  or  the  bottom  of  a 
conduction  band  near  its  Van  Hove 


singularity.  The  energy  of  the  top  of  the 
first  valence  band  (E^)  has  been 

chosen  to  be  zero  in  Table  1.  "Deg" 
refers  to  the  degeneracy  of  the  bands. 
Dm  (in  states/eV)  is  the  peak  value  of 
the  density  of  states  at  the  Van  Hove 
singularity;  it  was  used  to  calculate  the 
effective  mass  in* using  the  universal 
relation,  Eq.(7)  or  (9).  We  utilized  the 
effective  mass  of  charge  carriers  of 
SWCNT  (10,  0)  at  £v(1)  as  the  single 
reference  for  the  calculations  of  the 
effective  masses  (in*)  of  other 
SWCNTs.  As  a  comparison,  we  also 
calculated  the  effective  mass  m*  using 
the  results  of  the  conventional  method 
for  which  we  utilized  the  ab-initio 
electronic  structure  calculations.  We 
performed  self-consistent,  ab-initio 
density  functional  calculations  to  obtain 
the  electronic  structure  of  the  SWCNTs, 
using  a  linear  combination  of  atomic 
orbital  (LCAO)  method.610  Particularly, 
the  ab-initio  calculations  for  the 
electronic  structure  of  nonsymmorphic 
SWCNTs  (8,  4)  and  (10,  5)  are  very 
difficult  tasks,  since  there  are  112  and 
140  atoms  in  their  tubule  unit  cells, 
respectively.  The  calculated  results  of 
the  effective  masses  in* of  the  charge 
carriers  in  SWCNTs  utilizing  the 
universal  relation  Eq.  (7)  or  (9)  agree 
very  well  with  those  from  the 
conventional  method  of  ab-initio 
electronic  structure  calculations,  as 
demonstrated  in  Table  1. 
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Table  1.  The  calculated  effective  masses  of  charge  carriers  in  various  1-D  SWCNTs, 
where  m*  is  the  calculated  effective  mass  from  ab-initio  density  functional 
computations;  in*  is  obtained  using  the  universal  relation  Eq.  (7)  or  (9);  and  Dm  is  the 
peak  value  of  the  density  of  states  near  Van  Hove  singularity. 


Band  Em  Deg 

m*  (ab-initio) 

Dm 

in* (Eq.  7  &  9) 

(eV) 

(m0) 

(states/eV) 

(m0) 

SWCNT  (10,  0) 

0.0  2 

0.10 

0.73 

0.10 

0.95  2 

0.10 

0.72 

0.10 

SWCNT  (13,  0) 

£v(1) * III IV  0.0  2 

0.09 

0.66 

0.08 

£C(I)  0.75  2 

0.09 

0.66 

0.08 

SWCNT  (17,  0) 

£v(1>  0.0  2 

0.11 

0.73 

0.10 

£c0)  0.54  2 

0.11 

0.72 

0.10 

SWCNT  (22,  0) 

Ev0)  0.0  2 

0.057 

0.51 

0.05 

£c(1)  0.44  2 

0.056 

0.51 

0.05 

SWCNT  (8,  4) 

E™  0.0 

0.13 

0.76 

0.11 

£<U)  0.96 

0.13 

0.76 

0.11 

SWCNT  (10,  5) 

£<U)  0.0 

0.11 

0.79 

0.12 

£°'2)  0.74 

0.11 

0.79 

0.12 
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Abstract 

We  report  calculated,  electronic  properties  of  wurtzite  zinc  oxide  (w-ZnO).  Unlike  many 
previous  theoretical  works,  our  linear  combination  of  atomic  orbital  (LCAO) 
calculations,  implemented  following  the  Bagayoko,  Zhao,  and  Williams  (BZW)  method, 
employed  a  local  density  approximation  (LDA)  potential  to  obtain  band  gaps  and  other 
results  in  agreement  with  experiments.  We  discuss  the  band  structures,  the  calculated 
direct  band  gap  (3.2  to  3.3  eV),  and  the  electron  effective  mass  (0.25  to  0.26  mo). 

PACS  Numbers  71.20.Nr,  71.15.Ap,  71.15.Mb 

I.  Introduction  and  Motivations 

Nowotny  et  al  [1]  and  Thomas  [2]  provided  some  of  the  earliest  experimental  data  on 
wurtzite  zinc  oxide  (w-ZnO),  the  former  on  the  lattice  parameters  and  the  latter  on  the 
band  gap.  A  representative  set  of  measurements  of  the  band  gap  of  bulk  w-ZnO  shows  a 
robust  agreement  on  its  low  temperature  value  of  3.44  ±  0.06  eV.  [1-7]  This  range,  i.e., 
±  0.06  eV,  is  much  larger  than  experimental  uncertainties;  it  simply  reflects  differences 
that  a  traceable  to  those  of  actual  lattice  constants,  measurement  techniques  and 
temperatures,  and  other  factors.  Some  of  these  factors  (impurities,  defects)  are  associated 
with  the  growth  conditions  of  the  samples.  Srikant  and  Clarke  [8]  appear  to  have  settled 
the  debate  on  the  room  temperature  band  gap  of  w-ZnO.  They  utilized  different 
measurement  techniques  to  arrive  at  a  value  of  3.30  eV.  These  techniques  included 
reflection  and  transmission  absorption,  ellipsometric  spectroscopy,  Fourier  transform 
infrared  (FTIR)  spectroscopy,  and  photoluminescence.  This  multitude  of  measurement 
techniques  avoided  limitations  reported  to  afflict  some  previous  investigations  that 
reported  a  room  temperature  band  gap  of  bulk  w-ZnO  0.1  and  0.2  eV  below  the  currently 
accepted  value  of  3.30  eV.  Other  groups  reported  photoemission  [9-11]  results  for  the 
band  structure,  the  electrons  effective  mass  [12]  and  the  bulk  modulus  [13]. 

Several  groups  [14-20]  have  reported  measured  band  gaps  for  films  of  w-ZnO  grown  by  a 
variety  of  techniques.  These  growth  methods  included  spray  pyrolisis  [14],  pulsed  laser 
deposition  [16,18],  cathodic  electrodeposi.tion  in  aqueous  solutions  [17],  ion  layer  gas 
reaction  (ILGAR)  and  rf  magnetron  sputtering  [19],  and  plasma  enhanced  chemical  vapor 
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deposition  (PECVO).  [20]  Li  and  Coworkers  [20]  discussed  the  variation  of  the  lattice 
parameters  of  films  with  changes  in  growth  temperatures.  They  also  noted  that  despites 
these  variations,  the  lattice  parameters  (a  and  c)  of  the  wurtzite  films  are  within  3.4x10 3 
Angstroms  from  the  corresponding  values  for  the  bulk  material.  The  reported  band  gaps 
.for  the  films  are  within  the  range  of  -0.14  to  +0.06  from  3.44  eV.  Studenikin  et  al  [14] 
found  a  gap  of  3.30  eV,  after  annealing  at  400  C,  for  all  their  films  grown  by  spray 
pyrolisis.  This  result  was  confirmed  by  Pauporte  and  Lincot  [17]  for  films  grown  by 
catholic  electrodeposition  in  aqueous  solutions.  The  band  gap  of  their  films  as  deposited 
was  3.5  eV.  These  authors  [17]  noted  that  growth  conditions  affect  the  gap.  In 
particular,  the  presence  of  nitrate  or  perchlorate  ions  is  reported  to  lead  to  a  band  gap  of 
3.30  eV.  Using  both  transmission  and  photoluminescence  techniques,  Muth  et  al. 
obtained  band  gaps  of  3.40  eV  and  3.45  eV  at  295  K  and  77  K,  respectively.  Within  the 
understandable  differences  due  to  growth  stress,  thermal  expansion  mismatch,  and 
impurities  [8],  the  above  band  gaps  for  w-ZnO  films  are  comparable  to  those  for  high 
quality  bulk  w-ZnO  samples. 

In  contrast  to  the  above  agreement  between  experiments,  theoretical  calculations 
reported  values  of  the  band  gap  that  cover  a  wide  spectrum  and  mostly  disagree  with 
measurements.  The  Hartree  Fock  (HF)  calculations  [21-24]  led  to  band  gaps  above  11 
eV,  more  than  three  times  the  measured  value.  When  some  correlation  effects  are  taken 
into  account,  the  modified  HF  calculations  of  Jaffe  et  al  [21]  produced  a  gap  of  8.61  eV. 
Empirical  pseudopotential  [25]  and  empirical  KKR  [26]  calculations  predictably  led  to 
3.5  and  3.30  eV,  respectively.  These  values  are  close  to  the  experimental  gap  of  3.44  eV. 
The  fitting  involved  in  these  empirical  approaches  can  adjust  some  parameters  so  as  to 
reproduce  the  experimental  band  gap.  Local  density  approximations  (LDA)  potentials, 
depending  on  the  computational  methods  [23,  27-33],  led  to  values  between  0.23  eV  and 
2.26  eV  for  the  band  gap  of  w-ZnO.  Table  1,  further  below,  identifies  some  of  the 
various  methods  and  applicable  lattice  parameters.  The  computational  approaches  include 
the  pseudopotential  (PP),  the  augmented  plane  wave  (APW),  the  full-potential  linearized 
augmented  plane  wave  (FLAPW),  and  the  linearized  muffin-tin  orbital  (LMTO)  methods. 
The  Green  function  and  screened  coulomb  approximations  (GWA)  calculations  did  not 
totally  resolve  the  above  failure  of  ab-initio  calculations  to  obtain  the  measured  band  gap 
of  w-ZnO.  The  GWA  work  of  Oshikiri  and  Aryasetiawan  [31],  for  two  different  sets  of 
lattice  constants,  obtained  4.28  (4.06)  and  3.45  (3.63)  eV  respectively,  with  the  values 
between  parentheses  resulting  from  using  additional  higher  energy  product  basis  for 
improved  accuracy.  The  GWA  calculations  of  Usuda  et  al  [30]  obtained  a  gap  of  2.44  eV 
with  both  the  LMTO  and  LAPW  methods.  This  value  is  1  eV  below  the  measured  value 
of  3.44  eV. 

The  gross  disagreement  between  ab-initio  calculations  and  experiment,  as  far  as  the  band 
gap  is  concerned,  basically  precludes  an  extensive  discussion  of  calculated  electron 
effective  masses,  inasmuch  as  the  latter  are  determined  by  the  curvature  of  the  band  in  the 
immediate  vicinity  of  the  lowest  energy  in  the  conduction  band.  Indeed,  difficulties 
responsible  for  the  disagreement  on  the  gap  are  also  expected  to  affect  not  just  the 
location  of  the  lowest  energy  in  the  conduction  band  but  also  the  curvature  of  said  band. 
With  the  above  understanding,  we  note  that  Oshikiri  et  al  [34]  reported  quasiparticle  (and 
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LDA)  effective  masses  of  0.25  (0.21),  0.25  (0.21),  and  0.21  (0.18)  mo  in  the  kx,  ky,  and  k2 
directions  respectively,  where  the  values  between  parentheses  resulted  from  their  LDA 
calculations.  These  results  were  obtained  for  experimental  lattice  constants  of  3.24961  A 
,  5.20653  A  and  u  parameter  of  0.345.  Lambrecht  et  al.  [33]  noted  that  their  calculated 
values  of  0.23  mo  and  0.21  mo,  in  the  parallel  and  perpendicular  directions,  are  likely  to 
be  underestimates.  Xu  and  Ching  [35]  reported  0.37,  0.28,  and  0.32  mo  in  the  T-K,  T-A, 
and  T-M  directions  in  the  Brillouin  zone,  for  lattice  constants  given  in  Table  1.  These 
authors  utilized  the  orthogonalized  LCAO  formalism  and  employed  a  local  density 
potential  plus  additional  exchange  derived  with  Wigner’s  interpolation.  The  1966 
measurements  by  Baer  [12]  found  the  electron  effective  mass  for  w-ZnO  to  be  0.24 
±  0.02  mo.  The  above  LDA  results  for  the  electron  effective  mass  are  uniformly  lower 
than  this  experimental  value. 

The  above  overview  of  previous  theoretical  results,  with  the  general  disagreement 
between  calculations,  on  the  one  hand,  and  between  calculations  and  the  firmly 
established  experimental  data,  on  the  other  hand,  provides  the  motivation  for  this  work. 
The  wide  spread  of  the  LDA  results  for  the  band  gap  of  w-ZnO,  from  0.23  to  2.26  eV, 
suggests  that  limitations  other  than  those  genuinely  attributable  to  the  potential  affected 
these  results.  Since  1998,  Bagayoko,  Zhao,  and  Williams  have  introduced  a 
computational  method  [36-39]  that  clearly  showed  that  LDA  potentials  are  no  longer 
known  to  lead  to  woeful  underestimates  of  the  band  gap  of  semiconductors.  Indeed, 
using  the  Rayleigh  theorem,  these  authors  introduced  a  new  form  of  convergence  for 
variational  computations  that  utilize  a  basis  set.  The  need  to  seek  methodically  [36-37]  a 
basis  set  whose  size  is  converged  [i.e.,  the  optimal  basis  set]  vis  a  vis  the  description  of 
the  occupied  states,  they  argued,  straightforwardly  arises  from  (a)  the  use  of  only  the 
wave  functions  of  the  occupied  states  to  construct  the  charge  density,  the  potential,  and 
the  Hamiltonian  in  the  iterative  process  and  (b)  the  Rayleigh  theorem  that  asserts  the 
possible,  continuing  lowering  of  some  unoccupied  energies  or  bands  when  basis  sets 
larger  than  the  optimal  one  are  utilized.  This  extra-lowering  of  some  unoccupied 
energies,  even  though  the  physics  (i.e.,  the  Hamiltonian)  is  no  longer  changing,  was 
identified  as  a  basis  set  and  variational  effect  that  has  afflicted  most  of  the  previous 
calculations,  including  those  utilizing  the  generalized  gradient  approximation  (GGA)  of 
density  functional  theory  (DFT)  [40-41]  and  the  GWA  approximation  [42].  Hence,  the 
clear  motivation  of  this  work  is  to  utilize  properly,  as  per  the  BZW  method,  a  local 
density  potential  to  describe  the  electronic  and  related  properties  of  w-ZnO.  In  the 
remaining  of  this  article,  we  briefly  present  our  method  and  discuss  our  results,  which  are 
basically  in  agreement  with  experiment,  particularly  as  compared  to  previous,  theoretical 
ones. 


II.  Method 

Our  computational  method  has  been  extensively  described  in  previous  publications  [36- 
39].  It  is  characterized  by  the  use  of  a  local  density  approximation  (LDA)  potential  [43] 
as  parameterized  by  Vosko  et  al  [44],  the  well-known  linear  combination  of  Gaussian 
orbitals  (LOGO),  and  the  rigorous  application  of  the  BZW  method  in  carrying  out  the 
calculations.  We  utilized  a  program  package  developed  and  refined  over  decades  [45-46], 
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Our  calculations  are  non-relativistic  and  are  performed  at  zero  temperature.  They 
commence  with  the  LCAO,  self-consistent  calculations  of  the  electronic  energy  levels  of 
the  atomic  or  ionic  species  that  are  present  in  the  system  under  study.  For  wurtzite  ZnO, 
these  species  are  Zn  (Zn2+)  and  O  (O2').  Specifically,  the  calculations  for  the  neutral 
atoms  are  utilized  in  _preliminary  investigations  of  the  solid  state.  The  output  of  the 
preliminary  calculations  is  employed  to  compute  ionic  charges,  if  any,  in  the  solid.  We 
found  zinc  and  oxygen  elements,  in  the  solid,  to  be  closer  to  Zn2+  and  O2'  than  to  neutral 
Zn  and  O  or  Zn+1  and  O1'.  Consequently,  ab-initio  calculations  were  performed  for  Zn2+ 
and  O2'  in  order  to  obtain  the  “atomic”  basis  sets  to  be  employed  in  the  actual  solid  state 
calculations. 

These  self-consistent,  solid  state  calculations,  as  per  the  BZW  method,  started  with  the 
minimal  basis  (MB)  set,  i.e.,  the  one  just  enough  to  account  for  all  the  electrons  in  the 
system.  Calculation  II  utilized  a  basis  set  comprised  of  the  MB  set  plus  orbital(s) 
representing  the  next  excited  level(s)  in  the  ionic  (or  atomic)  species  in  the  system.  The 
occupied  energies  from  Calculations  I  and  II  were  compared  numerically  and  graphically. 
As  is  generally  the  case,  they  were  different.  Hence,  Calculation  HI  was  performed  with  a 
basis  set  including  that  of  Calculation  II  as  augmented  with  orbitals  representing  the  next 
excited  levels  of  the  ionic  species  in  the  system,  i.e.,  Zn2+  and  O2".  As  discussed  below, 
the  occupied  energies  from  Calculations  II  and  III  are  the  same,  within  computational 
uncertainties.  Hence,  the  basis  set  of  Calculation  II  is  our  optimal  basis  set,  i.e.,  the  one 
utilized  for  all  results  reported  here.  If  the  Rayleigh  theorem  did  not  apply  and  if  the 
wave  functions  of  the  occupied  states  were  not  the  only  ones  utilized  in  the  interactive 
process  to  construct  the  charge  density  and  the  Hamiltonian,  then  the  results  of 
Calculation  HI,  with  sinking  unoccupied  bands,  would  have  been  the  physical  answer  on 
the  account  of  basis  set  completeness  considerations.  The  actual  size  of  the  basis  set 
naturally  depends  on  the  types  of  functions,  with  plane  wave  approaches  entailing  many 
more  basis  functions  than  those  utilizing  exponential  or  Gaussian  functions. 

Computational  details  germane  to  the  replication  of  our  work  follow.  ZnO  (or  zincite) 
possesses  a  hexagonal  lattice  in  the  space  group  C4<5V.  There  are  four  atoms  per  unit  cell 
in  the  (2b)47  positions  as  follows:  Zn:  (0,  0,  0),  (1/3,  2/3,  1/2);  O:  (0,  0,  u),  (1/3,  2/3,  1/2  + 
u).  The  self-consistent  computations  were  performed  for  different  sets  of  lattice 
parameters  as  shown  in  Table  1.  For  a  given  calculation,  we  employed  a  mesh  of  24  k- 
points,  with  proper  weights,  in  the  irreducible  Brillouin  zone.  Our  criterion  for  self- 
consistency  rested  on  the  convergence  of  the  potential  to  a  difference  around  10‘5 
between  two  consecutive  iterations.  Approximately  60  iterations  were  needed  to  reach 
self-consistency.  The  computational  error  for  the  valence  charge  was  about  0.0002  for  52 
electrons.  For  the  implementation  of  the  BZW  method,  our  first  self-consistent 
calculation  utilized  the  atomic  orbitals  of  Zn2+(ls2s3s4s2p3p3d)  and  02‘(ls2s2p). 
Calculation  II  included  the  (4p)  orbital  on  zinc  while  Calculation  III  further  added  the  3p 
orbital  on  Oxygen.  The  comparison  of  the  occupied  states  from  Calculation  II  and  III 
showed  that  the  convergence  of  the  size  of  the  basis  set,  vis  a  vis  the  description  of  the 
occupied  states,  was  reached  with  Calculation  II.  The  basis  set  of  this  calculation,  the 
optimal  basis  set,  was  employed,  for  the  results  discussed  below. 
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III.  Results 


Some  of  our  calculated  band  structures  are  shown  in  Figures  1  and  2.  Figure  1  illustrates 
the  dramatic  differences  between  our  results  and  those  or  previous  LDA  and  other 
calculations.  Essentially,  there  are  no  major  differences  between  our  calculated  occupied 
bands  and  some  previously  reported  ones.  For  the  unoccupied  bands,  however,  our  results 
are  distinctively  far  those  of  previous  calculations.  In  Figure  1,  the  lowest  unoccupied 
energy  at  the  T  point  is  shifted  downward  by  approximately  2.2  eV  as  compared  to  its 
value  obtained  with  the  optimal  basis  set.  Clearly,  the  differences  between  the  low 
energy  conduction  bands  from  calculation  II  and  III  are  far  from  those  expected  from  a 
rigid  shift.  Hence,  they  cannot  be  corrected  with  simplistic  "scissors"  operations. 

The  above  differences  between  our  calculated  conduction  bands  and  those  in  previous 
reports  explain  the  reason  our  calculated  band  gaps  of  3.2  and  3.278  eV,  at  lattice 
constants  specified  in  Table  1  and  Figures  1  and  2,  are  very  close  to  experimentally 
measured  values.  The  differences  between  our  results  from  the  calculation  at  the 
experimental  lattice  constants  in  Figure  2  and  the  one  just  below  it,  in  Table  1,  are  mainly 
due  to  the  difference  in  the  value  of  u,  given  that  the  lattice  parameters  are  not  far  apart. 
Schroer  and  coworkers  found  a  similar  situation  in  their  LMTO  and  GWA  calculations. 
[27]  The  value  of  u,  understandably,  strongly  influences  that  of  the  band  gap. 

Our  calculated,  total  (DOS)  and  partial  density  of  states  (pDOS)  in  Figures  3  and  4 
naturally  reflect  the  differences  discussed  above  for  the  electronic  energy  bands.  In 
particular,  the  location  of  the  lowest  conduction  bands  with  respect  to  the  Fermi  energy  in 
these  figures  is  drastically  different  from  that  in  most  of  the  previous  theoretical 
investigations.  The  inset  in  Figure  3  illustrates  our  concept  of  a  practical  band  gap.  [37] 
Essentially,  due  to  instrumental  sensitivity,  analysis  techniques,  and  related  uncertainties, 
an  experimental  work  could  find  values  of  the  band  gap  from  the  theoretical  minimum  of 
3.278  to  close  to  4.0  eV.  This  possibility,  practically  non-existent  in  highly  sensitive 
photoluminescence  measurements,  increases  to  a  likelihood  if  simple  optical  absorption 
is  the  sole  measurement  technique.  The  determination  of  the  onset  of  the  band  edge 
absorption  is  a  source  of  large  uncertainties,  in  addition  to  the  fitting  often  involved  in  the 
determination  of  the  band  gap. 

The  electron  effective  masses  from  our  calculations  exhibit  a  slight  anisotropy,  as 
expected  in  a  wurtzite  structure.  The  effective  mass  is  a  measure  of  the  curvature  of  the 
calculated  bands.  The  agreement  between  calculated  and  measured  effective  masses 
indicates  an  accurate  determination  of  the  shape  of  the  bands.  We  calculated  the 
effective  masses  of  the  n-type  carriers  of  ZnO,  using  the  electronic  structure  from 
calculation  II  (the  solid  line  in  Fig.  1).  Near  the  bottom  of  the  lowest  conduction  band  at 
the  gamma  point,  we  obtained  0.254,  0.260,  and  0.264  mo  in  the  T-A,  T-K,  and  T-M 
directions,  respectively,  for  the  lattice  constants  in  Figure  1.  For  the  conduction  band  in 
Figure  2,  the  corresponding  values  are  0.257,  0.258,  and  0.257  mo  respectively.  These 
values  agree  very  well  with  the  measured  value  of  0.24  mo,  within  the  experimental 
uncertainty  [12]  of  0.02  mo,  and  with  the  experimental  value  of  0.275  m0.  [48]  The  above 
excellent  agreement  with  experiment  and  the  underestimates  from  previous  LDA 
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calculations  are  graphically  illustrated  in  Figure  1.  Indeed,  the  drastic  extra-lowering  of 
the  low  lying  conduction  band  (dotted  lines  in  Fig.  1),  in  calculations  employing  basis 
sets  larger  than  the  optimal  one,  results  in  a  significant  decrease  of  the  effective  mass.  In 
contrast,  the  lowest  lying  conduction  band  at  T,  as  obtained  with  the  optimal  basis  set 
(solid  lines  in  Fig  1),  leads  to  larger  values  of  the  electron  effective  mass. 

IV.  Conclusion 

We  have  performed  ab-initio,  self-consistent  calculations  of  the  electronic  energy  bands, 
the  density  of  states  (DOS),  and  of  the  electron  effective  masses  of  wurtzite  zinc  oxide 
(w-ZnO).  The  implementation  of  the  Bagayoko,  Zhao,  and  Williams  (BZW)  procedure, 
within  the  linear  combination  of  atomic  orbitals,  led  to  occupied  and  unoccupied  energy 
bands  respectively  similar  and  drastically  different  from  results  of  previous  calculations. 
Our  calculated  band  gap  of  3.278  eV,  at  an  experimental  lattice  constant,  is  only  5% 
smaller  than  the  measured  value  of  3.44  eV.  This  work  strongly  indicates  that  the 
previous  LDA  band  gaps  that  are  34%  to  93%  smaller  than  the  experimental  value  suffer 
from  the  basis  set  and  variational  effect  described  above  and  possibly  from  other 
computational  limitations;  hence,  their  gross  disagreement  with  experiment  should  no 
longer  be  ascribed  to  a  failure  of  LDA.  Our  calculated  electron  effective  masses  are  in 
excellent  agreement  with  experiment,  within  experimental  uncertainties. 
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Table  1.  Calculated,  fundamental  band  gaps  (Eg,  in  eV)  of  w-ZnO,  along  with 
pertinent  lattice  constants  in  Angstroms,  compared  to  experiment.  Numbers  in 
the  last  four  columns,  for  a  given  row,  are  from  the  reference  cited  in  that  row. 


Local  Density 
Approximation 
(LDA) 
Potentials 


Computational 

Method 


LCAO-BZW 
(Present  work) 


Zn2+  Pseudopotential 
Zn12+  Pseudopotential 


Pseudopotential 


APW 


LAPW 


FLAPW 


LMTO 


LMTO 


LMTO 


LMTO  (cubic,  a=4.57  A) 


FP-LMTO 


a  (A) 


3.2501 


3.2496 


3.2530 


3.2700 


3.25 


3.25 


3.23 


3.253 


3.2427 


3.2496 


3.253 


c(A) 


5.2071 


5.206 


5.2130 


5.1800 


5.21 


5.21 


5.18 


5.2129 


5.1948 


5.2065 


5.2129 


u 


0.3817 


0.345 


0.3817 


0.381 


0.3817 


0.3826 


0.345 


0.3825 


0.382 


Eg  (eV) 


3.20 


3.28 


3.21 


3.22 


2.13 


0.23 


0  23b 


1 .40 


0.77° 


0.93 


TTiT 


0.97' 


0.78° 


2.269 


1.813 


TT 


LDA  &  Additional 
Exchange _ 


Orthogonalized  LCAO 


3.249 


5.207 


0.345 


0.88' 


LDA+SIC 


Pseudopotential 


3.29 


5.29 


3.77 


Xa  Potential 


Pseudopotential 


3.25 


5.21 


1.58 


GWA 


LMTO 


3.2427 


5.1948 


0.3826 


4.28  and 
4.06a 


GWA 


LMTO 


3.2496 


5.2065 


0.345 


3.45  and 
3.63a 


GWA 


LMTO 


3.253 


5.2129 


0.3825 


2.44 


GWA 


LAPW 


3.253 


5.2129 


0.3817 


2.44d 


Model  GWA 


FLAPW 


4.23 


Empirical  Pseudopotential  (EMP) 


3.5J 


Empirical  KKR 


3.30 


Experiment  (Bulk,  Low  Temperatures) 


3.44 


Experiment  (Bulk,  Room  Temperature) 


3.30 


Reference  27 
dReference  30 
Reference  32 
'Reference  25 


“Reference  28 
Reference  23 
“Reference  33 
Reference  26 


Reference  29 


Reference  31 
'Reference  35 
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Figure  1.  Electronic  band  structures  of  wurtzite  zinc  oxide  (w-ZnO)  from  Calculations  II 
and  III  of  the  BZW  method.  The  solid  lines  show  the  results  of  calculation  II  and  the 
dashed  lines  represent  the  bands  from  calculation  III.  Calculation  II  gives  the  final 
electronic  structure  of  ZnO  with  the  optimal  basis  set.  The  Fermi  levels  from  the  two 

o 

calculations  are  superimposed.  The  lattice  constants  are  a-  3.2501  Angstroms  (A), 
c  =  5.2071  A,  and  u  =  0.3817 .  The  band  gap  is  3.2  electron  volts  (eV). 
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Figure  2.  Calculated  band  structure  of  wurtzite  zinc  oxide  (w-ZnO),  as  obtained  with 

o  o 

BZW  method,  for  experimental  lattice  constants  of  a-  3.2496  A,  c-  5.206  A,  and 
u  =  0.345 .  The  band  gap  is  3.278  eV. 
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DOS  (states/eV  cell) 


Figure  3.  Calculated  total  density  of  state  (DOS)  for  w-ZnO,  as  derived  from  the  bands 
shown  in  Fig.  2.  The  inset  indicates  that  experimental  apparatus  or  techniques  that  are 
not  very  sensitive  could  find  band  gaps  much  larger  than  the  theoretical  minimum  of 
3.278  eV. 


Energy  (eV) 
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Figure  4.  Calculated  partial  density  of  states  (pDOS)  for  w-ZnO  at  the  experimental 
lattice  constants  of  a-  3.2496  A,  c  =  5.206  A,  and  u  =  0.345  .  The  graph  illustrates  the 
influence  of  the  zinc  d  and  the  oxygen  p  in  the  upper  valence  bands.  Zinc  s  dominates  in 
the  lowest  conduction  band. 
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1.  INTRODUCTION  AND  MOTIVATIONS 

Several  review  papers1-4  have  discussed  the  properties 
and  applications  of  wurtzite  and  zinc-blende  indium  nitride 
(c-InN).  In  particular,  the  value  of  the  band  gap  of  wurtzite 
InN  (w-InN)  has  recently  attracted  much  interest3-5  due  to 
seemingly  conflicting  findings  from  experimental  investiga¬ 
tions.  The  current  and  potential  applications  of  InN  based 
semiconductor  devices  certainly  warrant  a  rapid  resolution  of 
the  unsettled  issues  relative  to  these  important  materials.1-5 
Indeed,  light  emitting  diodes  (LEDs),  laser  diodes  (LDs),  and 
photodiodes  (PDs),  over  a  wide  range  of  energy,  from  ultra¬ 
violet  to  infrared  wavelengths,  are  some  of  these  applications 
of  InN  based  semiconductors.  The  direct  band  gaps  of  high 
quality  wurtzite  InN  films  are  reported3-5  to  be  from  0.65  to 
1  eV,  depending  on  carrier  concentration  and  other  sample 
characteristics.  The  band  gap  of  cubic  InN  (c-InN)  is  ex¬ 
pected  to  be  in  this  range  or  slightly  below  it.  Hence,  InN 
films  could  be  critical  for  the  fabrication  of  high  speed  LDs 
arid  PDs  for  optical  communication  system's.3  Unlike  w 
-InN,  the  most  stable  phase  of  InN  in  ordinary  conditions, 
much  remains  to  the  known  about  c-InN.  This  situation  is 
partly  due  to  the  serious  difficulties  associated  with  the 
growth  of 

c-InN.3,4 

These  difficulties  are  apparent  in  the  work  of  Yamamoto 
et  al.6  who  grew  c-InN  on  GaAs  and  a-Al203  substrates. 
They  reported  the  appearance  of  w-InN  when  the  film  thick¬ 
ness  was  over  0.05  /zm,  and  at  thicknesses  over  0.2  /zm,-  the 
c-InN  films  grown  on  GaAs  were  completely  covered  by 
hexagonal  indium  nitride.  The  films  grown  on  sapphire  con¬ 
tained  columnar,  fibrous  structures.  Unlike  Yamamoto  et  al.6 
who  utilized  metalorganic  vapor-phase  epitaxy  (MOVPE) 
Tabata  et  al.,1  Chandrasekhar  et  al. ,8  and  Cimalla  et  al.9  em¬ 
ployed  molecular  beam  epitaxy  (MBE)  to  grow  c-InN  films. 
These  authors  reported  values  of  the  lattice  constant  of 
c-InN  of  4.97±  0.01,  4.980,  and  4.986  A,  respectively.  We 
are  unaware  of  reports  of  experimental  investigations  of  the 
electronic  and  elastic  properties  of  c-InN.  Hence,  an  aim  of 
this  work  is  to  predict  electronic  and  related  properties  of  c 
-InN,  including  the  band  gap. 

The  need  for  this  work  is  partly  underscored  by  the  dis¬ 
agreement  between  previous  theoretical  findings.10-22  Previ¬ 
ous  local  density  approximation  (LDA)  calculations  that  em¬ 

0021-8979/2004/96(8)/4297/5/$22.00 


ployed  the  pseudopotential  (PP)  method  produced  negative 
band  gaps  ranging10,11,13-13,19  from  -0.18  to  -0.40  eV.  Other 
LDA  calculations,11,12,16,19  using  variations  of  the  linearized 
augmented  plane  wave  (LAPW), 11,12  linear  muffin-tin  orbital 
(LMTO),16  and  the  atomic  sphere  approximation  (ASA) 
(Ref.  19)  obtained  c-InN  band  gaps  of  0.08-0.48  eV.  The 
generalized  gradient  approximation  (GGA),  ’  within  the 
pseudopotential  approach,  led  to  a  gap  value13  of  -0.55  eV. 
Self-interaction  corrections  to  LDA,  quasiparticle  (QP)  ap¬ 
proaches,  and  exact  exchange  calculations  reported  c-InN 
band  gaps  of  0.43-1.40  eV  as  shown  in  Table  I  below.  The 
very  recent  empirical  pseudopotential  result  of  Fritsch  et 
al.12  is  0.592  eV  for  the  band  gap  of  c-InN.  Our  focus  here  is 
mainly  on  the  LDA  results  that  are  very  small  or  negative  for 
the  band  gap  of  c-InN.  This  situation  constitutes  another  mo¬ 
tivation  for  this  work,  besides  the  lack  of  experimental  data 
on  the  electronic  and  related  properties  of  c-InN.  Our  LDA 
calculations  recently  resolved  the  controversy  that  was  sur¬ 
rounding  the  band  gap  of  wurtzite  InN.23  In  particular,  our 
LDA  calculations,  within  the  Bagayoko,  Zhao,  and  Williams 
(BZW)  implementation  of  the  linear  combination  of  atomic 
orbitals  (LCAO)  method,  obtained  a  w-InN  band  gap  of 
0.88  eV,  in  very  good  agreement  with  measurements  from  a 
recent  series  of  experimental  investigations.23 

Further,  utilizing  features  of  the  calculated  density  of 
states  (DOS),  we  showed"  the  possibility  of  obtaining  a  gap 
as  large  as  2  eV  if  optical  absorption  is  the  only  measure¬ 
ment  technique  utilized.  Difficulties  in  precisely  determining 
the  band  edge,  analysis  techniques,  and  related  uncertainties 
explain  this  assertion.  Our  work  on  w-InN  added  to  a  series 
of  articles  according  to  which  it  is  no  longer  correct  to  state 
that  local  density  approximation  woefully  underestimates  the 
band  gaps  of  semiconductors.  Hence,  and  in  the  light  of  the 
mostly  negative  values  of  the  band  gap  of  c-InN,  there  exists 
a  compelling  rationale  for  an  implementation  of  the  LDA 
that  does  not  suffer  from  an  effect24-26  believed  to  be  mostly 
responsible  of  the  dismal  underestimation  of  band  gaps.  In¬ 
cidentally,  the  band  gap  problem  was  actually  a  symptom  of 
a  more  general  and  unrecognized  problem  stemming  from 
unoccupied  levels  or  bands  that  are  affected  by  the  effect 
referenced  above  and  recalled  below. 

©  2004  American  Institute  of  Physics 
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TABLE  I.  Experimental  and  theoretical  lattice  constants  (a,  in  angstroms)  for  c-InN,  along  with  the  calculated  values  of  the  bulk  modulus  gigapascal  and  of 
the  fundamental  band  gap  (eV).  Results  in  the  last  three  columns,  for  a  given  row,  are  from  the  reference  cited  in  that  row. 

Computational  Method  a(k)  B(GPa)  £g(eV ) 


LCAO-BZW  (Present  work) 

5.017 

4.98 

145 

+0.65 

+0.74 

4.95 

145 

-0.36“ 

5.004 

140 

-0.40b 

Pseudopotential  Method  (PP) 

4.95 

4.97d 

145" 

4.932 

140 

-0.35' 

Local  Density  Approximation  (LDA)  Potentials 

4.788 

155f 

-0.18s 

LAPW 

4.94 

145c 

Full  potential  LAPW 

5.03 

138 

-0.11h 

-0.48h 

Full  Potential  LMTO 

4.92 

139* 

-0.4' 

Atomic  sphere  approximation  (ASA) 

-0.11 

+0.02  and  +0.08 

Generalized  gradient  approximation  (GGA) 

PP 

5.06 

5.109 

120c 

118 

-0.55b 

LDA  plus  self-interaction  correction  (SIC) 

5.05“ 

+0.43“ 

QP  Calculation 

PP 

+0.52“ 

QP+SIC 

PP 

+1.31“ 

DFT  Exact  Exchange 

+1.46 

DFT,  SX 

ASA 

+1.31 

Estimate  of  the  bulk  modulus  of  zinc-blende  indium  nitride  (c-InN)  using  elastic  properties  of  wurtzite  InN 

1371 

Empirical  Pseudopotential  Calculations  (EMP) 
Experimental:  Measured  lattice  constants 

4.97±0.01“ 

+0.592"* 

4.98° 

4.986p 

“Reference  10.  -* 

‘Reference  16. 

"Reference  13. 

^Reference  17. 

"Reference  11. 

"Reference  19. 

"Reference  15. 

‘Reference  20. 

“Reference  18. 

“"Reference  22. 

Reference  21. 

"Reference  7. 

Reference  14. 

"Reference  8. 

"Reference  12. 

^Reference  9. 

II.  METHOD  AND  COMPUTATIONAL  DETAILS 

We  performed  zero  temperature,  nonrelativistic  calcula¬ 
tions  of  the  electronic  and  related  properties  of  c-InN.  Our 
ab  inito ,  self-consistent  calculations  employed  the  local  den¬ 
sity  approximation  potential  of  Ceperley  and  Alder27  as  pa¬ 
rametrized  by  Perdew  and  Zunger.28  As  stated  above,  we 
used  the  LCAO.  The  feature  distinguishing  our  computa¬ 
tional  method  from  the  previous  investigations  noted  above 
consists  of  our  implementation  of  the  BZW  procedure.  In  so 
doing,  we  started  the  calculations  for  c-InN  with  a  minimal 
basis  set.  We  subsequently  performed  several  other  self- 
consistent  calculations  with  larger  and  larger  basis  sets.  The 
basis  set  for  any  of  these  calculations  was  obtained  by  aug¬ 
menting  the  one  for  the  previous  calculation  with  the  orbital 
describing  the  next  excited  level  of  the  atomic  or  ionic  spe¬ 
cies  present  in  the  system.  The  occupied  bands  of  a  given 
calculation  are  compared  to  those  of  the  previous.  These 
comparisons,  particularly  for  the  first  two  calculations,  often 
show  differences  (in  numerical  values,  branching,  or  curva¬ 
ture).  This  process  continues  until  the  occupied  energies 


from  a  calculation  are  equal,  within  computational  uncertain¬ 
ties,  to  their  corresponding  ones  from  the  calculation  that 
follows  it.  Then,  the  output  of  the  former  calculation  pro¬ 
vides  the  physical  description  of  the  material  under  study  and 
the  related  basis  set  is  dubbed  the  optimal  basis  set.  Accord¬ 
ing  to  the  Rayleigh  theorem,26  some  of  the  unoccupied  bands 
from  the  latter  may  be  lower  than  their  counterpart  from  the 
former.24-26'29,30 

The  above  additional  lowering  is  the  basis  set  and  varia¬ 
tional  effect  inherently  associated  with  variational  calcula¬ 
tions  of  the  Rayleigh-Ritz  type.  In  the  iterative  process,  the 
use  of  the  wave  functions  for  the  occupied  states  only  in  the 
construction  of  the  charge  density,  and  hence  the  potential 
and  the  Hamiltonian,  ensures  the  exhaustion  of  the  account¬ 
ing  for  the  physical'interactions  when  the  occupied  energies 
converge  vis-a-vis  the  size  of  the  basis  set.  As  fully  ex- 
plained  elsewhere,  "  however,  some  unoccupied  ener¬ 
gies  will  continue  to  be  lowered  as  the  basis  set  is  increased 
beyond  the  optimal  one.  The  sizes  of  the  minimal  and  opti¬ 
mal  basis  sets  vary  vastly  with  the  type  of  functions  in  the 
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TABLE  II.  Eigenvalues  (eV),  along  high  symmetry  points,  for  zinc-blende 
indium  nitride  (c-InN)  as  obtained  from  LDA-BZW  calculations  for  a 
=5.017  A,  the  theoretical  equilibrium  value.  The  Fermi  energy  of 
-0.216  87  eV  is  set  to  zero  in  the  table. 


L 

r 

X  r 

K 

-15.5192 

-14.8427 

-15.6493 

-15.5949 

-14.7559 

-14.8427 

-14.7301 

-14.7034 

-14.7559  ' 

-14.8427 

-14.5254  ' 

-14.5837 

-14.4740 

-14.5336 

-14.5254 

-14.5739 

-14.4740 

-14.5336 

-14.4463 

-14.4865 

-11.7624 

-14.1793 

-11.3085 

-11.2078 

-5.5713 

0.0000 

-4.7495 

-4.5900 

-0.8335 

0.0000 

-2.2093 

-3.3631 

-0.8335 

0.0000 

-2.2093 

-1.7848 

4.0316 

0.6536 

4.1816 

5.6110 

8.4339 

9.7615 

6.8978 

7.1493 

10.6892 

9.7615 

11.8922 

10.5553 

10.6892 

9.7615 

11.8922 

11.6130 

atomic  orbitals  (i.e.,  exponential,  Gaussian,  plane  wave  func¬ 
tions)  and  with  other  features  of  the  calculations  (i.e., 
pseudopotential,  LAPW,  etc).  As  described  above,  the  BZW 
method  is  applicable  in  most  computations  that  utilize  the 
LCAO  formalism.  The  LCAO  program  package  we  employed 
in  this  work  was  developed  many  years31-33  before  the  intro¬ 
duction  of  the  BZW  method  in  1998.  We  provide  below 
computational  details  germane  to  a  replication  of  our  calcu¬ 
lations. 

Zinc-blende  InN  is  a  member  of  the  III-V  family.  The 
atomic  wave  functions  of  the  ionic  states  of  In1+  and  N1- 
were  obtained  from  self-consistent  ab  initio  calculations.  The 
radial  parts  of  the  atomic  wave  functions  were  expanded  in 
terms  of  Gaussian  functions.  A  set  of  even-tempered  Gauss¬ 
ian  exponents  was  employed  for  In  with  a  minimum  of 
0.1400  and  a  maximum  of  0.2300  X  106,  with  16  Gaussian 
functions  for  the  s'  and  p  states  and  14  for  the  d  states  of  In1+. 
The  two  largest  exponents  were  not  included  in  the  descrip¬ 
tion  of  the  d  state.  Similarly,  a  set  of  even-tempered  Gauss¬ 
ian  exponents  was  utilized  to  describe  N1-,  with  a  minimum 
of  0.1242  and  a  maximum  of  0.1365  X  105.  Both  the  s  and  p 


TABLE  III.  Eigenvalues  (eV),  along  high  symmetry  points,  for  zinc-blende 
indium  nitride  (c-InN),  as  obtained  from  LDA-BZW  calculation  for  a 
=4.98  A.  The  Fermi  energy  -0.216  58  eV  is  set  to  zero  the  table. 


L 

r 

X 

K 

-15.5977 

-14.8586 

-15.7269 

-15.6720 

-14.7642 

-14.8586 

-14.7370 

-14.7084 

-14.7642 

-14.8586 

-14.5160 

-14.5781 

-14.4583 

-14.5212 

-14.5160 

-14.5685 

-14.4583 

-14.5212 

-14.4279 

-14.4714 

-11.7856 

-14.3473 

-11.3252 

-11.2093 

-5.7063 

0.0000  • 

-4.8350 

-4.6807 

-0.8675 

0.0000 

-2.2907 

-3.4730 

-0.8675 

0.0000 

-2.2907 

-1.8515 

4.1748 

0.7383 

4.2140 

5.6870 

8.5382 

9.9148 

7.1020 

7.3237 

10.8738 

9.9148 

12.1018 

10:7467 

10.8738 

9.9148 

12.1018 

11.8177 

FIG.  1.  Calculated  LDA-BZW  band  structure  of  zinc-blende  Indium  Nitride 
(c-InN)  at  the  theoretical  equilibrium  lattice  constant  of  5.017  A,  as  ob¬ 
tained  with  BZW  optima]  basis  set.  The  Fermi  level  (-0.216  87  eV)  is  set  to 
zero  in  the  figure. 

functions  were  expanded  in  terms  of  13  Gaussian  orbitals. 
The  minimal  basis  set  comprised  atomic  orbitals  representing 
In1+  (\s2s2p3s3p3d4s4p4dSs5p)  and  N1-  (Is2s2p).  A  mesh 
of  28 k  points  in  the  irreducible  Brillouin  zone,  with  proper 
weights,  was  used  in  the  self-consistent  iterations.  The  com¬ 
putational  error  for  the  valence  charge  was  about 
-0.001  449  79  for  36  electrons.  The  self-consistent  poten¬ 
tials  converged  to  a  difference  around  10-5  after  about  60 
iterations. 

III.  RESULTS  AND  DISCUSSIONS 

Table  I  shows  our  calculated,  theoretical  equilibrium 
band  gap  of  0.65  eV,  at  a  lattice  constant  of  5.017  A,  and  the 
value  of  0.74  eV  at  an  experimental  lattice  constant  of 
4.98  A.  Tables  II  and  III  contain  the  calculated  energies,  at 
some  high  symmetry  points  in  the  Brillouin  zone,  for  the  two 
lattice  constants  given  above.  Figures  1-3  exhibit  the  energy 
bands,  for  the  theoretical  equilibrium  lattice  constant 
(5.017  A),  and  the  related  total  (DOS)  and  partial  (pDOS) 
densities  of  states.  The  curve  of  the  total  energy  versus  the 
lattice  constant  is  shown  in  Fig.  4. 

The  lack  of  experimental  data,  except  for  the  lattice  con¬ 
stant,  precludes  an  extensive  discussion  of  these  results.  Our 
predicted  equilibrium  lattice  constant  is  within  0.6%  from 
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FIG.  2.  Total  DOS  for  zinc-blende  indium  nitride  (c-InN)  as  obtained  with 
the  bands  shown  in  Fig.  1.  The  inset  illustrates  our  definition  of  a  “practical 
band  gap.” 
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FIG.  3.  pDOS  for  zinc-blende  indium  nitride  (c-InN)  as  obtained  with  the 
bands  shown  in  Fig.  1.  The  dominance  of  nitrogen  p  at  the  top  of  the 
valence  band  is  obvious  in  this  graph. 
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the  latest  experimental  one  of  4.986  A.  Our  calculated  bulk 
modulus  of  145  GPa  is  close  to  the  results  from  some  other 
LDA  calculations10,11  and  disagrees  with  findings  around 
120  GPa  from  GGA  calculation11,13  as  shown  in  Table  I.  The 
tables  for  the  calculated  energies  are  expected  to  provide 
useful  comparisons  for  future  experimental  investigations. 
Our  calculated  LDA-BZW  band  structure  is  drastically  dif¬ 
ferent  from  the  findings  of  most  of  the  previous 
calculations — as  far  as  the  unoccupied  bands  are  concerned. 
The  large  differences  between  our  calculated  band  gaps  and 
previous  LDA  results  in  Table  I  are  a  direct  consequence  of 
the  differences  between  the  respective  unoccupied  bands. 

The  empirical  pseudopotential  (EMP)  result22  of 
0.592  eV  is  relatively  close  to  our  findings  of  0.65  and 
0.74  eV.  We  cannot  draw  much  fundamental  significance 
from  this  fact,  however,  given  that  this  EMP  result  was  ob¬ 
tained  using  model  potential  parameters  whose  derivation 
entailed  fitting  to  data  that  included  a  band  gap  of  0.59  eV. 
For  the  purposes  of  application,  however,  this  closeness  por¬ 
tends  much  importance.  Indeed,  we  expect  potential  param¬ 
eters  derived  from  fitting  to  our  data  to  lend  themselves  to 
credible  and  practical  descriptions  of  electronic,  optical,  elas¬ 
tic,  and  structural  properties  of  materials.  This  assertion  is 
partly  supported  by  the  versatility  and  relative  ease  of  em¬ 
pirical  pseudopotential  calculations. 

As  in  the  case  of  w- InN, 23  our  results  do  not  show  any 
indication  of  an  overestimation  of  the  p-d  repulsion  by  LDA 
potentials.  This  overestimation  was  believed10  to  be  the 
source  of  the  very  small  or  negative  band  gaps  by  pushing 
the  top  of  the  valence  band,  dominated  by  p  states,  to  higher 
energies.  According  to  our  findings,  it  is  rather  the  extra  low¬ 
ering  of  the'  bottom  of  the  conduction  band  that  produces 
LDA  band  gaps  that  are  negative  or  very  small23,25  if  LCAO 
type  computations  do  not  utilize  the  BZW  approach  to  avoid 
it  while  ensuring  the  adequacy  (or  convergence)  of  the  basis 
set  for  the  description'  of  occupied  states.  In  fact,  the  basic 
derivation  of  the  ground  state  theory  that  is  the  original  den¬ 
sity  functional  theory'34,35  implicitly  suggests  such  an  ap¬ 
proach,  notwithstanding  the  need  to  account  for  the  redistri¬ 
bution  of  electrons  in  molecular  to  solid  environments  with 

36 

judicious  polarization  and  diffuse  orbitals. 
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FIG.  4.  Total  energy  of  (c-InN)  vs  the  lattice  constant.  The  total  energy  at 
the  equilibrium  lattice  constant  of  5.017  A  is  1949.4215  eV. 

As  per  their  definition,  effective  masses  provide  a  mea¬ 
sure  of  the  quality  (i.e.,  curvature)  of  band  structures.  Our 
LDA-BZW  calculations  found  electron  effective  masses,  at 
the  bottom  of  the  conduction  band,  of  0.065m0,  0.066m0,  and 
0.06 6m0  in  the  T-L,  F-X,  and  F-K  directions,  respectively, 
for  the  equilibrium  lattice  constant.  For  a  lattice  constant  of 
4.98  A,  the  corresponding  effective  masses  for  the  electrons 
are  0.07 6m0,  0.073mo,  and  0.073 m0  in  the  T-L,  F-X,  and 
T-K  directions,  respectively.  Our  calculated,  equilibrium, 
electron  effective  masses  are  very  close  to  the  0.066m0  from 
the  EMP  calculations  of  Fritsch,  Schmidt,  and  Grundmann.22 
This  agreement  supports  our  comment  above  relative  to  the 
potential  use  of  EMP  calculations  when  the  potential  param¬ 
eters  are  derived  in  part  by  fitting  to  LDA-BZW  results. 

IV.  CONCLUSION 

In  lieu  of  a  conclusion,  we  contend  that  our  ab  initio, 
self-consistent  LDA-BZW  calculations  have  predicted  elec¬ 
tronic,  structural,  and  elastic  properties  of  cubic  InN  (c 
-InN)  in  the  zinc-blende  structure.  It  is  hoped  that  experi¬ 
mental  investigations  will  follow  in  the  near  future.  It 
emerges  from  this  work  that  theoretical  efforts  should  be 
directed  to  the  determination  of  actual  limitations  of  LDA 
and  of  other  approaches  as  opposed  to  echoing  the  chorus 
now  known  to  be  untenable,  as  per  tire  physical  interactions, 
and  that  ascribes  to  LDA  a  woeful  underestimation  of  band 
gaps  that  is,  as  per  the  Rayleigh  theorem,  straightforwardly  a 
consequence  of  a  basis  set  and  variational  effect. 
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We  report  the  calculated  band  gap  of  wurtzite  indium  nitride.  Our  ab  initio  computations  employed 
a  local-density  approximation  (LDA)  potential  and  the  linear  combination  of  Gaussian1  orbital 
formalism.  The  implementation  of  the  ab  initio  Bagayoko,  Zhao,  and  Williams  method  [Phys.  Rev. 
B  60,  1563  (1999)]  led  to  a  LDA  band  gap  of  0.88  eV,  in  excellent  agreement  with  recent 
experiments.  We  also  present  calculated  density  of  states  (DOS)  and  the  electron  effective  mass  at 
the  bottom  of  the  conduction  band.  Our  DOS  curves  indicate  that  an  experiment  could  find  values 
of  the  band  gap  up  to  2  eV,  depending  on  the  sensitivity  of  the  apparatus,  the  interpretation  of 
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BACKGROUND  AND  MOTIVATIONS 

The  1998  review  articles  of  Orton  and  Foxon1  and  of 
Ambacher  extensively  described  the  technological  impor¬ 
tance  and  potential  of  group  Ill-nitrides,  including  InN.  The 
ternary  alloy  In^xGaxN  offers  a  wide  range  of  band  gaps 
that  depend  on  X.  Until  recently,  the  lower  limit  of  this  range 
was  believed1,2  to  be  about  1.8-2.0  eV.  This  limit  was  the 
previously  measured  band  gaps  of  wurtzite  InN  as  obtained 
from  the  absorption  studies  of  mostly  polycrystalline  films 
grown  by  dc  discharge,3  various  sputtering  techniques,4-8  or 
metal-organic  vapor-phase  epitaxy  (MOVPE).9  The  above 
measured  band  gaps  of  wurtzite  InN,  i.e.,  1. 8-2.0  eV,  are  in 
disagreement  with  findings  of  recent  experiments  discussed 
below. 

The  2000  study  of  the  band  bowing  of  InAIN  by  Yama- 
gushi  et  al.10  strongly  indicated  that  the  band  gap  (Eg)  of 
wurtzite  InN  should  be  less  than  1.9  eV.  Subsequent  mea¬ 
surements  of  the  band  gap  of  InN  films,  mostly  grown  by 
molecular-beam  epitaxy  (MBE),  agreed  on  a  range  of 
0.7- 1.0  eV.  The  high  quality  of  the  samples  in  these  experi¬ 
ments  was  partly  utilized  to  explain  the  large  difference  be¬ 
tween  the  early3-9  and  latter10-  7  groups  of  experiments.  In- 
ishima  et  alu  found  the  value  of  Eg  to  be  0.89  and  1.46  eV 
for  InN  films  at  respective  electron  concentrations  of  5 
X  1019  and  2  X  1020  cm-3.  The  multipronged  study  by  Wu  et 
al ,12  determined  the  band  gap  of  w-InN  by  optical  absorp¬ 
tion,  photoluminescence,  and  photomodulated  reflectance 
techniques.  The  room-temperature  band  gap  was  found  to  be 
between  0.7  and  0.8  eV  for  their  high-quality  films  that  were 
grown  by  molecular-beam  epitaxy.  These  authors  further  es¬ 
tablished  that  no  noticeable  feature  could  be  found  in  the 
absorption  and  photomodulated  reflectance  spectra  to  indi¬ 
cate  a  band  gap  between  1 .9  and  2.0  eV.  The  band-gap  bow¬ 
ing  of  InGaN  also  pointed  to  a  value13  of  Eg  near  0.8  eV  for 
wurtzite  InN.  Even  though  Matsuoka  et  al. 13  studied  samples 
of  wurtzite  InN  (w-InN)  grown  by  metal-organic  vapor- 
phase  epitaxy,  one  of  the  growth  methods  of  the  early  group 
of  experiments,  they  found  that  Eg  for  their  w-InN  samples 
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was  between  0.7  and  1.0  eV.  The  2004  work  of  Arnaudov  et 
al.  on  InN  films  grown  by  MBE  or  MOVPE  gives  a  value 
of  0.692±0.002  eV  for  the  low-temperature  (2  K)  band  gap 
of  high-quality  samples  with  relatively  low  electron  concen¬ 
trations  (i.e.,  7.7  X  1017-6x  1018  cm-3).  These  authors  em¬ 
ployed  near-band-edge  photoluminescence.  The  above  agree¬ 
ment  of  recent  experiments  on  Eg  values  between  0.7  and 
1.0  eV,  depending  on  temperature,  pressure,  and  carrier  con¬ 
centration,  is  in  stark  contrast  with  the  large  range  of  the 
following  theoretical  findings. 

Early  empirical  pseudopotential  model  calculations18,19 
(EPM)  basically  reproduced  the  gap  of  2.0  eV  reported  by 
the  above  first  group  of  experiments.  The  fitting  involved  in 
the  method  practically  permits  the  reproduction  of  any  se¬ 
lected  feature  of  the  electronic  structure.  The  2003  EPM 
work  of  Fritsch  et  al.  not  only  reproduced  an  earlier  experi¬ 
mental  band  gap  of  2.1 1  eV,  but  also  reported  values  of  2.58 
and  2.59  eV,  for  wurtzite  InN,  from  their  anisotropic  and 
isotropic  pseudopotential  calculations. 

Local-density  approximation  (LDA)  results  for  the  band 
gap  of  wurtzite  InN,  including  those  of  authors  who  also 
performed  quasiparticle-  ’  and  other  calculations,  range  ' 
from  -0.4  to  0.43  eV.  The  pseudopotential22’24-26  (PP)  and 
the  full  potential  linearized  augmented  plane-wave21,23 
(FLAPW)  calculations  mostly  reported  negative  values  of 
the  band  gap,  from  -0.4  to  -0.19  eV.  The  full  potential  lin¬ 
ear  muffin-tin  orbital23  (FP-LMTO)  and  atomic  sphere  ap¬ 
proximation  (ASA)  approaches  led  to  values  between  0.2 
and  0.43  eV. 

A  pseudopotential  calculation  within  the  generalized  gra¬ 
dient  approximation25  (GGA)  produced  a  band  gap  of 
-0.37  eV.  Several  authors  reported  the  band  gap  of  wurtzite 
InN  as  obtained  by  quasiparticle22,25  (QP)  or  other  computa¬ 
tions  that  employed  various  schemes  that  either  go  beyond 
the  local-density  approximation  ’  ’  or  alter  LDA 
potentials.  ’  ’  In  particular,  the  pseudofunction"  ap¬ 
proach,  self-interaction  correction  (SIC),22  and  self- 
interaction  with  relaxation  corrections  (SIRC),24  respectively, 
led  to  band  gaps  of  1.3,  0.58,  and  1.6  eV.  The  GW  method 
led  to  results22,27  of  0.74  and  0.5  eV.  The  GW-SI  produced  a 
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sap  of  1.5  eV.  The  linear  combination  of  Gaussian  orbital 

^  OQ 

(LCGO)  calculations  of  Xu  and  Ching  employed  a  LDA 
potential  plus  additional  exchange  effects  obtained  with 
Wigner  interpolation  formula.  They  reported  a  band  gap  of 
■  1.02  eV.  The  LDA+C  work  of  Wei  et  al.21"  augmented  the 
LDA  potential  with  atom-dependent  corrections  and  obtained 
a  band  gap  of  0.85  eV  with  an  uncertainty  of  0.1  eV. 

Clearly,  the  above  state  of  theoretical'  calculations  is  am¬ 
ply  sufficient  to  motivate  the  present  work.  A  recent  experi¬ 
mental  development30  exacerbated  the  urgency  of  our  theo¬ 
retical  investigations.  Indeed,  the  concordant  findings  of  the 
second  group  of  experiments,  described  above,  have  just 
been  questioned  by  Shubina  et  al.30  in  2004.  These  authors 
showed  the  possibility  of  explaining  the  spectra  of  this  sec¬ 
ond  group  of  experiments  not  in  terms  of  near-band-edge 
absorption  or  emission,  but  rather  as  manifestations  of 
surface/gap  states  or  Mie  resonances  due  to  In  precipitates  in 
InN.  They  noted  the  poor  thermal  stability  of  InN  and  the 
low  In  vapor  pressure  to  support  their  position.  Specifically, 
for  high-quality  InN  films  grown  at  T>  500  °C  by  MBE  or 
MOVPE,  they  utilized  x-ray  diffraction  to  observe  tetragonal 
In  precipitates  in  hexagonal  InN.  Their  very  low-temperature 
(0.35  K)  measurements  by  thermally  detected  optical  ab¬ 
sorption  (TDOA)  indicate  the  possibility  of  having  a  band 
gap  around  1.4  eV.  The  authors  refrained  from  making  a 
definitive  statement  as  to  the  location  of  the  fundamental 
band  gap,  despite  their  observations,  due  to  the  difficulties  in 
separating  the  intraband  'absorption  in  InN  from  Mie  reso¬ 
nances. 

METHOD  AND  COMPUTATIONAL  DETAILS 

Our  calculations  employed  the  Ceperley  and  Alder31 
local-density-functional  potential  as  parametrized  by  Perdew 
and  Zunger32  in  an  ab  initio  LCGO  formalism.  We  also  em- 
ployed  the  GGA  potential  of  Perdew  and  co-workers  which 
goes  beyond  LDA.  The  calculations  are  nonrelativistic  and 
for  zero  temperature.  A  key  distinction  between  this  work 
and  previous  calculations  resides  in  our  utilization  of  the  ab 
initio  Bagayoko,  Zhao,  and  Williams34'35  (BZW)  method  in 
carrying  out  the  self-consistent  computations.  Introduced34  in 

oo 

1998  and  extensively  employed  thereafter,  this  method 
rigorously  avoids  a  basis  set  and  variational  effect  inherently 
associated  with  all  variational  calculations  that  utilize  both 
the  linear  combination  of  atomic  orbital  (LCAO)  and  a 
charge  density  obtained  only  with  the  wave  functions  of  the 
occupied  states.  The  effect,  a  consequence  of  the  Rayleigh 
theorem,  consists  of  a  lowering  of  unoccupied  energy  levels 
or  bands,  which  is  not  ascribable  to  a  physical  interaction, 
after  the  total  convergence  of  the  charge  density,  the  poten¬ 
tial,  and  of  the  occupied  energy  levels  with  respect  to  the  size 
of  the  basis  set.34,30 

BZW  self-consistent  calculations  begin  with  the  mini¬ 
mum  basis  set,  i.e.,  one  that  is  just  large  enough  to  include  all 
the  orbitals  for  the  occupied  levels  in  the  atomic  species  that 
are  present  in  the  system  under  study.  This  initial  basis  set,  in 
our  case,  is  derived  from  calculations  for  the  atomic  or  ionic 
species  present  in  the  finite  or  infinite  system  under  investi¬ 
gation.  Completely  self-consistent  calculations  with  this  ba¬ 


sis  set  are  generally  not  expected  to  describe  the  system  cor¬ 
rectly.  Indeed,  they  generally  do  not  account  for  electron 
redistribution  in  the  system  as  opposed  to  that  in  neutral 
atoms  or  ionic  species.  A  second  self-consistent  calculation  is 
performed  with  a  basis  set  that  includes  the  minimum  basis 
set  as  augmented  with  the  orbital  representing  the  lowest 
excited  levels  in  the  atomic  (or  ionic  species).  The  occupied 
energy  levels  or  bands  of  calculations  I  andTl  are  compared. 
One  generally' finds  differences  between  these  energies  from 
the  two  calculations.  For  solids,  these  differences  are  in  nu¬ 
merical  values  and  the  branching  and  curvatures  of  the 
bands.  A  third  calculation  is  performed  by  augmenting  the 
basis  set  of  calculation  II  as  described  above.  If  a  comparison 
of  the  self-consistent,  occupied  energies  of  calculations  Et 
and  III  do  not  show  a  complete  agreement,  within  applicable 
errors,  the  process  of  augmenting  the  basis  set  is  continued. 
Eventually,  calculations  N  and  (N+l)  will  have  the  same 
occupied  energies.  Further,  the  charge  density  and  the  poten¬ 
tial,  from  these  two  calculations,  are,  respectively,  the 
same.34,35  Upon  the  convergence  of  the  variational  calcula¬ 
tions  with  respect  to  the  size  of  the  basis,  the  results  for 
calculation  N  are  selected  as  the  physical  ones.  One  gener¬ 
ally  finds,  indeed,  that  calculation  (N+l)  produced  some  un¬ 
occupied  energies  that  are  lower  than  the  corresponding  ones 
in  calculation  N.  These  differences  that  vary  vastly  with  the 
value  of  ic  in  reciprocal  space  are  manifestations  of  the  afore¬ 
mentioned  basis  set  and  variational  effect.34-38  The  basis  set 
for  calculation  N  is  called  the  optimal  basis  set. 

Other  ordinary  features  of  our  calculations  are  thor¬ 
oughly  described  in  the  literature  in  seminal  papers  that  re¬ 
ported'  on  the  development39-41  of  the  program  package  and 
its  recent  utilization.34'35,38  We  provide  below  computational 
details  germane  to  a  replication  of  our  ab  initio  work.  The 
experimental,3  low-temperature  wurtzite  lattice  parameters 
utilized  in  our  work  are  a=3.544A,  c=5.718  A,  and  u 
=  0.3790.  Preliminary  calculations  for  the  neutral  In  and  N 
atoms  and  the  solid  w-InN  indicated  that  the  self-consistent 
system  is  approximately  In0-s+N°'S“.  Consequently, 
we  performed  self-consistent  calculations  for  In1+  and  N1- 
to  obtain  the  input,  “atomic”  orbitals.  The  minimal  basis 
set  comprised  atomic  orbitals  representing  In1+ 
(\s2s2p3s3p3d4s4p4d5s5p)  and  N1-  (Is2s2p).  We  em¬ 
ployed  even-tempered  Gaussian  functions  in  the  construction 
of  the  atomic  orbitals.  The  s  and  p  orbitals  of  In  were  de¬ 
scribed  with  16  even-tempered  Gaussian  functions  with  the 
respective  minimum  and  maximum  exponents  of  0.14  and 
0.2300  X  106.  The  d  orbitals  for  In  entail  14  even-tempered 
Gaussian,  the  same  as  for  the  ones  s  and  p  at  the  exclusion  of 
the  two  largest  exponents.  Similarly,  the  s  and  p  orbitals  of  N 
were  constructed  using  13  even-tempered  Gaussian  expo-, 
nents  whose  smallest  and  largest  values  were  0.124  215  and 
0.136  000  5  X  105.  The  computations  were  generally  •  com¬ 
pletely  self-consistent  after  60  iterations.  The  convergence, 
for  a  given  self-consistent  calculation,  was  reached  when  the 
potential  did  not  change  by  more  than  10-5  from  one  itera¬ 
tion  to  the  next.  Upon  reaching  self-consistency,  the  error 
made  in  accounting  for  72  electrons  was  0.004  31. 
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FIG.  1.  LDA-BZW  band  structure  of  wurtzite  indium  nitride  (InN).  Lattice 
constants:  a=3.544  A,  c=5.718  A,  and  tt= 0.3790.  The  solid  lines  represent 
the  band  as  obtained  with  the  optimal  basis  set.  The  dotted  lines  show  an 
example  of  bands  resulting  from  calculations  with  basis  sets  larger  than  the 
optimal  one.  The  Fermi  level  is  set  to  zero.  The  occupied  energies  are  the 
same  in  the  two  calculations. 
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FIG.  3.  Partial  density  of  states  (pDOS)  for  wurtzite  InN  as  obtained  from 
the  bands  (solid  lines)  shown  in  Fig.  1.  The  dominance  of  nitrogen  p  at  the 
top  of  the  valence  band  is  obvious  in  this  graph. 
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RESULTS 

Our  results  for  the  energy  bands,  the  total  (DOS),  and 
the  partial  densities  of  states  (pDOS)  are,  respectively,  shown 
in  Figs.  1-3.  The  LDA-BZW  bands  are  very  similar  to  the 
ones  obtained  with  GGA-BZW.  The  salient  difference  be¬ 
tween  these  results  and  the  findings  of  previous  calculations 
resides  in  the  description  of  the  unoccupied  bands.  Indeed, 
these  bands,  as  clearly  explained  above  in  relation  to  the 
Rayleigh  theorem  and  the  BZW  method,  do  not  suffer  from 
the  extra  lowering  ascribed  to  the  basis  set  and  variational 
effect  elucidated  by  Bagayoko  and  co-workers.34~38  It  is  im¬ 
portant  to  note  that  the  referenced  lowering  of  some  unoccu¬ 
pied  bands  is  far  from  being  rigid,  i.e.,  it  is  highly  dependent 
on  the  k  point  and  is  greatest  for  free-electronlike  band  (para¬ 
bolic  or  not)  and  minimum  for  flat  ones  and  at  high  symme¬ 
try  points.  Table  I  provides  the  energies  at  high  symmetry 
points  in  the  Brillouin  zone.  It  is  hoped  that  these  data  will 
lend  themselves  to  further  comparison  of  our  results  with 
experimental  findings,  including  optical  transition  energies. 
The  LDA-BZW  calculated  band  gap  of  wurtzite  indium  ni- 
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FIG.  2.  Total  density  of  states  (DOS)  for  wurtzite  indium  nitride  (InN)  as 
obtained  with  the  bands  (solid  lines)  shown  in  Fig.  1.  The  inset  supports  our 
definition  of  a  “practical  band  gap.” 


tride  is  0.88  eV.  The  GGA-BZW  gap  is  0.81  eV.  These  val¬ 
ues  agree  very  well  with  the  second  and  recent  group  of 
experiments  described  above. 

Bagayoko  et  al ,42  predicted  the  LDA  band  gaps  of  0.65 
and  0.74  eV  for  InN  in  the  zinc-blende  structure  (c-InN),  for 
the  theoretical  equilibrium  (5.017  A),  and  a  low-temperature 
experimental  (4.98  A)  lattice  constants,  respectively.  These 
values  are  approximately  0.2  eV  smaller  than  the  above  find¬ 
ings  for  wurtzite  InN.  Since  the  publication  of  these 
predictions,42  we  have  come  across  the  recent  experimental 
findings43  that  somewhat  agree  with  them,  within  the  limita¬ 
tions  associated  with  a  high  density  of  defects  in  the  cubic 
layer,  hexagonal  inclusions  up  to  30%  to  the  surface,  and  a 
very  high  carrier  concentration  above  1019  cm-3.  Indeed,  Ci- 
malla  et  al.  reported,43  within  the  preceding  limitations,  the 
results  of  optical  investigations  that  indicated  an  absorption 
edge  below  0.6  eV.  The  authors43  recommended  confirma¬ 
tion  of  this  indication  with  a  high-quality  cubic  InN. 

Our  group34,35  introduced  the  concept  of  “practical  band 
gap.”  Essentially,  the  practical  band  gap  designates  a  range 
of  values  of  the  band  gap,  as  per  the  inset  of  the  total 
density-of-state  (DOS)  curve,  such  that  an  experimental  mea¬ 
surement  could  produce  results  anywhere  in  that  range,  de¬ 
pending  on  the  sensitivity  of  the  apparatus  and  techniques, 
the  fitting  and  interpretations  of  these  results,  and  related 
uncertainties.  Clearly,  the  inset  in  Fig.  2  shows  that  measure¬ 
ments,  particularly  optical  absorptions,  could  lead  to  values 
of  the  band-gap  ranging  from  the  theoretical  minimum  of 
0.88  to  2  eV.  We  certainly  contend  that  the  provided  defects, 
impurities,  and  other  factors  (i.e.,  electron  concentration)  do 
not  explain  the  relatively  large  band  gap  of  w-InN  as  found 
by  the  first  group  of  experiments,  then  the  above  concept  of 
practical  band  gap  does.  Unlike  the  three-pronged  approach 
of  Wu  et  a/.,12--4  using  highly  sensitive  techniques,  optical 
absorption  was  generally  the  lone  approach  utilized  in  the 
early  determinations  of  the  band  gap  of  w-InN.  Even  the 
2004  results  of  Shubina  et  al,  as  mentioned  earlier,  are 
subject  to  this  interpretation.  Indeed,  these  authors  identified 
a  kink  to  denote  the  position  of  the  band  edge  even  though 
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TABLE  I.  Eigenvalues  (eV),  at  some  high  symmetry  points,  for  wurtzite  indium  nitride  (w-InN)  as  obtained 
fromLDA-BZW  calculations  for  a=3.544  A,  c=5.718  A,  and  u=0.3790.  The  Fermi  energy  of -0.222  12  eV  is 
set  to  zero  in  the  table. 


A 

L 

M 

H 

K 

r 

-13.0972 

-11.5668  ' 

-11.2329 

-11.3021 

-11.2827 

-11.6988 

-3.25532' 

-5.25538 

-5.20231 

-4.89999 

-4.32501 

-5.61294 

-3.25532 

-5.25538 

-4.27059 

-4.89999 

-4.32501 

-0.89662 

-0.46069 

-1.62671 

-3.34893 

-3.54349 

-4.2695 

-0.89662 

-0.46069 

-1.62671 

-2.29475 

-3.34349 

-2.39761 

-0.0049 

-0.46069 

-1.53365 

-1.75651 

-1.23894 

-2.39761 

0 

-0.46069 

-1.53365 

-0.8656 

-1.23894 

.  -2.14019 

0 

3.366891 

4.570732 

4.748424 

6.125875 

6.740313 

0.878663 

3.367163 

4.570732 

6.175673 

6.125875 

6.740313 

3.662409 

the  absorption,  as  per  their  Fig.  2(c)  visibly  commenced 
around  0.9  eV  and  reached  a  plateau  at  2  eV,  in  accordance 
with  our  concept  of  practical  band  gap. 

The  LDA-BZW  effective  masses  for  the  electron  close  to 
the  bottom  of  the  lowest  conduction  band  at  gamma  show  a 
slight  anisotropy.  The  electron  effective  mass  is  around 
0.088 m0  in  the  plane  perpendicular  to  kz  (from  T  to  M  and  F 
to  K)  and  0.082mo,  in  the  direction  of  kz  (from  F  to  A),  in  the 
immediate  vicinity  of  the  gamma  (r)  point.  The  GGA-BZW 
calculations  led  to  values  of  the  electron  effective  mass  very 
similar  to  the  ones  for  LDA.  Specifically,  at  the  bottom  of  the 
lowest  conduction  band,  they  are  0.081,  0.077,  and  0.082me 
from  T  to  M,  T  to  K,  and  F  to  A  directions,  respectively. 
These  effective  masses  are  also  in  agreement  with 
experiment,14  taking  into  account  the  quality  and  other  fea¬ 
tures  of  the  samples  (i.e.,  temperature,  pressure,  electron 
concentrations).  Wu  et  a/.1-4  obtained  a  value  of  0.07  by  using 
a  band  gap  of  0.7  eV  in  their  analysis.  They  also  provided  a 
plot  of  the  electron  effective  mass  as  a  function  of  electron 
concentration. 

DISCUSSIONS 

For  a  discussion  of  our  results,  we  underscore  the  need 
for  specificity  relative  to  the  lattice  parameters,  sample 
quality30  and  growth  conditions,11  temperature, 12and  electron 
concentration11  when  comparing  results  from  different  ex¬ 
perimental  investigations.  In  particular,  the  significant  effect 
of  electron  concentration  on  the  band  gap11’44'46  could  ex¬ 
plain  the  0.75-0.8-eV  band  gap  found  by  Wu  et  al.12  whose 
samples  had  a  relatively  low  electron  concentration  of  5 
X  1018  cm-3.  Our  calculated  band  gap  of  0.88  eV  agrees  very 
well  with  the  experimental  result  of  0.89  eV  of  Inishima  et 
al.n  for  a  sample  with  an  electron  concentration  of  5 
X  1019  cm-3.  As  noted  above  in  connection  with  the  “practi¬ 
cal  band-gap,”  measurement  techniques,  the  analysis  of  data, 
and  related  uncertainties  are  critical.  Specifically,  we  reiter¬ 
ate  that  absorption  techniques  could  lead  to  a  value  of  the 
band  gap  of  w-InN  as  large  as  2  eV. 

Some  authors22  evoked  a  p-d  repulsion  to  explain  the 
small  (generally  negative)  band  gap  from  previous  LDA  cal¬ 
culations.  LDA  is  reported  to  overestimate  this  repulsion  and 
to  push  upward  the  top  of  the  valence  band.  Our  results 
provide  no  indication  of  such  a  limitation  of  LDA.  Rather, 


they  clearly  show  that  the  extra  lowering  of  the  bottom  of  the 
conduction  band,  on  account  of  the  Rayleigh  theorem, 
straightforwardly  leads  to  the  band-gap  problem  in  LDA  and 
other  variational  calculations  that  employ  a  basis  set.  This 
work  partly  responds  to  the  calls  by  Ghani  et  al.44  and  by 
Vurgaftman  and  Meyer45  for  revisiting  the  issue  of  the  band 
gap  of  w-InN  following  their  extensive  compilations  of  ex¬ 
perimental  and  theoretical  results  on  indium  nitride.  Ghani  et 
al.  discussed  the  formation  of  oxynitrides  (InON),  particu¬ 
larly  in  aged  or  annealed  polycrystalline  samples,  as  a  pos¬ 
sible  source  for  the  increase  of  the  band  gap  to  values  around 
2  eV.  These  polycrystalline  films  are  reported  to  be  suscep¬ 
tible  to  containing  a  high  density  of  oxygen  atoms  at  their 
grain  boundaries. 


CONCLUSION 

In  conclusion,  we  claim  to  have  resolved  the  controversy 
surrounding  the  band  gap  of  wurtzite  InN.  We  did  so  within 
the  local-density  approximation  (LDA)  with  the  BZW 
method.  Unlike  many  schemes  purporting  to  correct  limita¬ 
tions  of  local-density  approximation  and  that  come  with 
other  approximations  whose  effects  are  at  best  ill  understood, 
the  ab  initio  BZW  method  is  rigorously  anchored  on  a  math¬ 
ematical  theorem.  The  above  points  explain  the  reason  this 
communication  is  hoped  to  end  the  echoing  of  the  incorrect 
chorus  claiming  that  it  is  well  known  that  local-density  ap¬ 
proximation  woefully  underestimates  the  band  gap  of  semi¬ 
conductors  and  insulators.  In  the  above  sense,  and  even 
though  this  communication  is  devoted  to  wurtzite  InN,  the 
method  applies  to  other  calculations  that  employ  basis  sets  in 
a  variational  scheme  of  the  Rayleigh-Ritz  type,  as  is  the  case 
in  many  energy  level  (finite  systems)  or  band  calculations. 
An  utterly  important  point  associated  with  the  above  points 
stems  from  the  fact  that  the  proliferation  of  schemes  aimed  at 
“correcting”  limitations  of  LDA  has  started  to  look  like  that 
of  epicycles  in  -the  Ptolemaic  model  of  the  solar  system.  A 
seriously  unfortunate  consequence  of  this  situation  consists 
of  diverting  efforts  from  grasping  the  actual  capabilities  and 
limitations  of  LDA  and  those  of  schemes  aimed  at  going 
beyond  the  local-density  approximation  or  beyond  density- 
functional  theory  altogether. 
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Abstract 

In  this  work,  we  reexamined  some  of  the  mathematical  and  phys¬ 
ical  properties  of  the  ab-initio  LCAO  calculations  for  the  electronic 
structure  of  ZnSe,  Ge,  and  GaAs.  The  utilization  of  non-strongly  min¬ 
imal  systems  in  the  self-consistent  ab-initio  calculations  could  lead  to 
a  non-uniformity  in  approaching  the  solution  in  the  Ritz-process.  We 
performed  test  computations  for  the  electronic  structure  of  hydrogen 
atom.  We  have  proposed  that  an  optimum  basis  set  may  be  needed 
so  that  the  calculated  electron  density  is  converged  and  the  significant 
scattering  of  the  Ritz-coefficients  may  be  avoided.  We  have  applied 
the  new  method  to  the  calculations  of  the  electronic  structure  of  ZnSe, 

Ge,  and  GaAs.  Our  calculated  results  of  the  electronic  properties  agree 
well  with  experimental  data. 

1  Introduction 

Ab-initio  density  functional  calculation  is  a  very  powerful  tool  to  study  the 
physical  properties  of  materials,  including  semiconductors,  metals,  surfaces, 
interfaces,  and  others. [1,  2,  3,  4,  5,  6,  7]  However,  previous  density  functional 
calculations  in  local  density  approximation  (LDA)  for  semiconductors  or  in¬ 
sulators  often  underestimated  the  band-gaps  by  30-50  %.[8,  9,  10, 11]  Closely 
related  to  the  band-gap  problem,  the  calculated  effective  masses  of  electrons 
(n-type  carriers)  and  the  optical  properties  of  semiconductors  -from  previ¬ 
ous  LDA  computations-  also  disagree  with  experimental  results.  There  have 
been  some  theoretical  efforts  intended  to  address  these  problems,  including 
the  calculations  that  utilize  nonlocal,  energy-dependent,  non-Hermitian  self¬ 
energy  operators. [12,  13,  14,  15,  16,  17,  18]  Others  implemented  the  exact 
exchange  potentials  or  included  the  effects  of  core  states  via  the  exchange 
diagram. [19,  20,  21]  Aryasetiawan  and  Gunnarsson  reviewed  several  compu¬ 
tational  methods,  including  the  GW  method,  aimed  at  describing  excited 
state  properties. [17].  Johnson  and  Ashcroft  recently  utilized  simplified  ap¬ 
plications  of  the  GW  method  to  make  scissors-type  corrections  to  the  band 
gaps  of  semiconductors.  [18] 
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In  this  work,  we  reexamine  some  of  the  mathematical  properties  of  the 
self-consistent  ab-initio  density  functional  calculations  that  utilize  the  lin¬ 
ear  combination  of  atomic  orbital  (LCAO)  formalism.  We  discuss  the  Ritz- 
process  of  the  variational  calculation  that  was  implemented  in  the  LCAO 
calculations  and  the  related  mathematical  theorem.  The  identified  mathe¬ 
matical  properties  of  the  Ritz-process,  which  is  the  foundation  of  the  LCAO 
calculations,  may  be  indispensable  for  understanding  the  band-gap  and  re¬ 
lated  problems  for  semiconductors.  We  performed  test  computations  for  the 
electronic  structure  of  hydrogen  atom.  We  will  illustrate  some  of  the  mathe¬ 
matical  and  physical  properties  of  the  self-consistent  ab-initio  calculations  of 
the  Ritz-process  by  the  results  of  the  calculations  on  hydrogen  atom,  ZnSe, 
Ge,  and  GaAs. 

In  the  next  section,  we  summarize  the  computational  method  and  the  re¬ 
lated  mathematical  properties.  We  present  the  calculated  results  in  Section 
III,  followed  by  a  short  conclusion  in  Section  IV. 


2  Method 


2.1  The  Kohn-Sham  Equation  and  The  Ritz  Process  in  the 
LCAO  Calculations 


Our  calculations  of  the  electronic  structure  of  materials  are  based  on  the 
density  functional  theory  of  Hohenberg-Kohn  and  Kohn-Sham.  [1,  2]  The 
Schrodinger-like  equation,  which  is  also  known  as  the  Kohn-Sham  equation, 
of  many  electron  system  [1,  2,  3]  is 


r  ^  2 
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+  e2/T|^ri*'  +  nc}^-i  =  Bji*{j  (1) 


where  and  T ^  are  the  eigen-energy  and  eigen-function  of  the  i-th  elec¬ 
tronic  state  at  the  /c- point  in  the  Brillouin  zone;  fm  and  Zm  are  the  posi¬ 
tion  and  nuclear  charge  of  the  m-th  ion;  Ri  is  a  translational  vector;  Vxc 
is  the  exchange-correlation  potential  of  the  many-electron  system.  In  this 
work,  we  employed  a  non-local  density  functional  formalism  of  the  exchange- 
correlation  potentials  in  the  generalized  gradient  approximation  (GGA)  that 
was  developed  by  Perdew  and  Wang.  [22]  In  fact,  the  Kohn-Sham  equation 
is  a  non-linear  differential  equation,  since  the  electron  density  p{r)  depends 
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on  the  wave  functions  Tr.  of  the  occupied  electron  states,  i.  e., 

=  l*«l2  (2) 

where  is  the  Fermi  distribution  function.  The  solution  of  the  Kohn-Sham 
equation  has  to  satisfy  a  condition  that 

J  p(r)dv  =  N  (3) 

where  N  is  the  total  number  of  electrons  in  the  system.  In  this  work,  we  im¬ 
plemented  the  formalism  of  the  linear  combination  of  atomic  orbitals  method 
to  solve  the  Kohn-Sham  equation.  The  essentials  of  the  mathematical  prop¬ 
erties  of  the  calculations  are  those  of  the  Ritz-process  of  the  variational  cal¬ 
culations.  We  will  follow  the  notations  in  the  discussion  of  the  variational 
theory  according  to  S.  G.  Mikhlin[23]. 

In  the  LCAO  method, [24,  25,  26,  27,  28,  29,  30,  31,  32]  we  solve  the 
Kohn-Sham  equation  self-consistly.  We  expand  the  eigenfunction  Tr-  of  the 
effective  single  electron  state  of  the  many-atom  system  as  a  linear  combina¬ 
tion  of  the  atomic  wave  functions,  which  is  a  Ritz  process  as  termed  in  the 
variational  theory  according  to  Mikhlin’s  notations.  [23] 

=Y.C<*rn{ki)<j>«m{k,?)  (4) 

am 

where  Cam  are  the  expansion  coefficient  or  the  Ritz-coefficients;  (/>am{k,f) 
are  Bloch  wave  functions  which  are  expressed  in  terms  of  atomic  wave  func¬ 
tions  as 

4>am{k',iE)  —  ^lUam{r  —  Tm  —  R{)  (5) 

where  uam  is  the  atomic  wave  function  of  the  a-th  state  of  the  m-th  atom 
at  the  position  rTO  and  is  the  coordinate  element  in  the  energy  space  (or 
Hilbert  space).  The  solution  for  the  Ritz-coefficients  Cam  of  the  Kohn- 
Sham  equation  is  the  one  that  the  total  energy  functional  Et0tai  (p)  reaches 
its  minimum. 

2.2  The  Requirement  of  Strongly  Minimal  Systems  in  the 
Ritz  Process 

We  shall  study  the  countable  system  of  the  coordinate  elements  in  the 
Hilbert  space,  where 

{hra}  =  ^1)^2)  U3, ...,  Un,  ...  (6) 
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The  corresponding  Gramm  matrix  of  the  first  n-elements  of  { un }  is 


/  (tti,ui)  (ux,u2)  ...  (uuun)  \ 
{u2,ui)  {u2;u2)  ...  [u2,un) 

V(  un >  ■Ui)  ( u 

n?  u2)  ...  (u  n 5  'U-n  )/ 


(7) 


where  ( Ui,Uj )  can  be  an  overlap  matrix  element  in  the  LCAO  formalism. 
The  eigen- values  of  the  Gramm  matrix  can  be  written  in  increasing  order  as 

\ni<\n2<  ...  <  K  (8) 

The  system  (Eq.  (6))  is  called  strongly  minimal  in  the  energy  space  (Hilber 
space),  if 

m/A"  =  lirrin-^ooXi  >  0  (9) 

A  relevant  mathematical  theorem  for  the  Ritz-process  in  the  variational 
calculation  is  the  following:  In  order  for  a  Ritz-process  to  be  stable,  it  is 
necessary  and  sufficient  that  its  generating  coordinate  system  { Ui }  be  strongly 
minimal  in  the  corresponding  energy  space.  [23] 

In  the  LCAO  calculations,  the  coordinate  system,  which  is  the  basis  set  of 
the  atomic  orbitals,  is  not  strongly  minimal  as  required  in  the  mathematical 
theorem.  The  problem  is  more  obvious  in  some  cases  of  the  implementations 
of  the  basis  sets  such  as  6-311G,  6-31G,  or  5-31G.  In  those  cases,  some  of 
the  shell  functions  that  describe  the  atomic  states  are  not  orthogonalized 
even  on  a  single  atomic  site.  In  the  solid  state  calculations,  the  smallest 
eigen-value  of  the  overlap  matrix  could  be  negative  in  those  calculations. 
As  experienced  by  various  researchers,  the  numerical  instability  problems 
could  occur,  which  are  usually  referred  as  ” linear  dependency”  problems. 

In  our  implementations  of  the  LCAO  method,  the  atomic  wave  functions 
(or  atomic  orbitals)  are  obtained  from  ab-initio  atomic  computations.  The 
atomic  shell  functions  are  orthogonal  in  the  cases  of  free  atoms  or  ions. 
However,  the  basis  sets  (or  coordinate  systems)  of  the  atomic  orbitals  for  the 
solid  state  calculations  are  neither  orthogonal  nor  strongly  minimal  as  that 
required  by  the  above  mentioned  mathematical  theorem.  The  problem  is 
more  intricate  since  we  may  notice  the  computational  instabilities  in  the  case 
that  the  problem  becomes  very  severe,  but  we  may  not  readily  appreciate  the 
beginning  of  the  appearance  of  the  computational  instabilities.  Particularly, 
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the  onset  of  the  computational  instabilities  in  the  LCAO  method  may  not  be 
marked  by  an  occurrence  of  the  negative  eigen-value  of  the  overlap  matrix. 

Mikhlin  illustrated  the  problems  of  the  utilization  of  non-strongly  min¬ 
imal  systems  in  Ref.  [23]  by  several  examples.  Generally,  the  non-strongly 
minimal  systems  will  not  reach  a  solution  in  the  limit  of  n  — )■  oo  in  the 
Ritz-process.  For  finite  values  of  n,  some  of  the  non-strongly  minimal  sys¬ 
tems  may  approximate  the  Ritz-coefficients  with  a  tendency  to  stabilize 
non-uniformly  for  increasing  n.  Other  non-strongly  minimal  systems  may 
show  that  the  Ritz-coefficients  scatter  significantly.  The  coefficients  may 
change  abruptly  from  one  approximation  to  the  next.  One  may  not  know 
the  exact  solution  for  the  problem  and  has  the  difficulty  of  identifying  an 
estimate  for  the  error  in  the  approximate  solutions.  There  is  so  far  no  math¬ 
ematical  theorem  that  can  identify  exactly  the  properties  of  the  instabilities 
in  the  utilization  of  non-strongly  minimal  systems  in  the  Ritz-process.  The 
relevant  problems  may  include  those  such  as  the  condition  (s)  that  the  in¬ 
stabilities  may  occur,  the  behavior  of  the  instabilities,  and  the  method (s)  to 
identify  them. 

The  problems  became  more  complex  in  the  applications  of  the  numerical 
computations  of  the  electronic  structure  of  solid  state  materials,  since  most 
of  the  basis  sets  of  various  implementations  of  ab-initio  density  functional 
calculations  are  not  strongly  minimal.  The  true  plane-wave  basis  set  sat¬ 
isfies  the  condition  to  be  strongly  minimal.  However,  the  requirement  of 
the  computational  accuracy  and  the  number  of  the  true  plane- waves  in  the 
solid  state  calculations  are  usually  beyond  the  limitation  of  our  computer 
resource  at  present.  The  orthogonalized  plane  waves,  which  are  derived  by 
the  Schmitdt  process  to  be  orthogonal  to  the  core  wave  functions,  are  not 
strongly  minimal. 

One  may  propose  to  orthogonalize  the  atomic  orbitals  for  the  solid  state 
calculations.  However,  if  the  basis  set  is  not  strongly  minimal,  the  orthog- 
onalization  of  the  basis  set  in  solid  state  calculations  may  encounter  the 
similar  instability  problem.  The  orthogonalization  procedure  of  the  basis 
set  involves  the  computation  of  the  inverse  of  the  overlap  matrix.  If  the 
smallest  eigen-value  of  the  overlap  matrix  approaches  zero  or  negative,  an 
instability  problem  may  occur. 

2.3  Consideration  of  Some  Physical  Properties 

Although  the  true  plane-wave  basis  set  can  be  implemented  in  ab-initio 
density  functional  calculations  of  the  electronic  structure  of  solid  state  ma- 
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terials,  the  LCAO  method  is  a  much  more  efficient  method  in  solving  the 
Kohn-Sham  equation  of  many-atom  systems.  The  utilization  of  non-strongly 
minimal  systems  (or  basis  sets)  such  as  those  used  in  LCAO  calculations 
could  lead  to  either  a  divergent  solution  as  n  ->  oo,  or  could  approach  their 
solutions  non-uniformly.  The  non- uniformity  in  approaching  the  solutions 
of  the  Ritz-process  or  the  scattering  of  the  Ritz-coefficients  are  particularly 
severe  when  the  sizes  of  the  basis  sets  become  relatively  large.  This  can  be 
seen  in  the  examples  demonstrated  by  Mikhlin  and  were  observed  in  pre¬ 
vious  ab-initio  electronic  structure  calculations  utilizing  the  LCAO  method 
on  various  materials. [23,  30,  31] 

In  the  cases  that  the  ab-initio  LCAO  calculations  or  other  similar  ones 
are  needed,  we  may  utilize  the  considerations  of  some  physical  properties 
of  the  solid  state  materials  in  identifying  approximate  solutions  without 
serious  errors.  One  of  the  plausible  hypotheses  may  be  the  following:  If  the 
exact  solutions  of  the  atomic  wave  functions  for  the  atoms  in  the  solid  state 
are  known,  the  required  size  n  of  the  basis  set  is  then  determined  by  these 
atomic  wave  functions.  There  will  be  no  need  for  any  additional  atomic 
orbitals  (or  atomic  wave  functions)  in  the  corresponding  energy  space  of 
the  exact  solution.  However,  we  may  not  know  the  exact  solutions  of  the 
atomic  wave  functions  in  the  solid  state  in  the  practical  calculations  and  may 
have  to  use  an  augmented  basis  set  of  atomic  orbitals.  We  will  address  this 
problem  in  more  detail  in  the  next  section  when  we  present  the  calculation 
results.  The  difficulty  may  be  lessened  if  a  global  self-consistent  approach 
is  implemented.  We  first  calculate  the  atomic  wave  functions  of  free  atoms 
in  a  neutral  state.  We  use  these  atomic  wave  functions  as  the  basis  set  to 
start  the  ab-initio  calculation  of  the  electronic  structure  of  the  solid.  We 
then  utilize  the  electronic  structure  of  the  solid  state  to  estimate  the  charge 
transfer  among  the  constituent  atoms.  The  information  is  consequently 
utilized  to  identify  the  corresponding  ionic  state  for  a  further  refinement  of 
the  calculations.  The  global  self-consistent  procedure  continues  until  the 
solutions  converge. 

3  Results 

3.1  Test  Calculations  on  Hydrogen  Atom 

In  the  calculations  of  the  electronic  structure  of  ZnSe,  Ge,  and  GaAs,  we 
implemented  self-consistent  ab-initio  density  functional  computations  that 
utilize  the  linear  combination  of  atomic  orbital  formalism.  To  obtain  the 
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atomic  wave  functions,  we  performed  ab-initio  atomic  calculations  that  also 
utilized  density  functional  formalism.  The  radial  parts  of  the  atomic  wave 
functions  were  expanded  in  terms  of  a  linear  combination  of  Gaussian  or¬ 
bitals  (LOGO). 

The* electronic  energy  levels  and  related  wave  functions  of  hydrogen  atom 
are  exactly  known  from  the  analytical  solution  of  the  Schrodinger  wave  equa¬ 
tion.  We  can  test  the  numerical  performance  of  the  variational  computations 
of  the  Ritz-process  on  hydrogen  atom.  We  utilized  the  same  ab-initio  atomic 
computation  program  and  the  LOGO  formalism  to  calculate  the  electronic 
structure  of  hydrogen  atom.  The  Schrodinger  wave  equation  for  the  calcula¬ 
tion  of  the  electronic  structure  of  hydrogen  atom  is  a  one-electron  problem, 
instead  of  the  many-body  problems  such  as  those  for  the  studies  of  solid 
state  materials.  In  the  ab-initio  computation  of  the  hydrogen  atom  using 
the  LOGO  method,  we  did  not  include  the  density  functional  potentials 
in  the  Hamiltonian.  The  Schrodinger  wave  equation  for  hydrogen  atom  is 
then  a  linear  differential  equation,  instead  of  the  non-linear  differential  equa¬ 
tion  such  as  the  density  functional  Kohn-Sham  equation.  We  compared  the 
results  of  the  numerical  computations  of  hydrogen  atom  with  the  exact  solu¬ 
tion.  In  the  numerical  computations,  we  utilized  basis  sets  of  even-tempered 
Gaussian  exponentials  with  a  minimum  of  0.10  x  10“ 2  and  a  maximum  of 
0.15  x  104  in  atomic  unit.  We  used  various  sizes  of  the  basis  sets  of  the 
Gaussian  functions. 

In  Fig.  1,  we  present  the  calculated  electronic  energy  levels  of  hydrogen 
atom  versus  the  total  numbers  of  the  utilized  Gaussian  functions  (Gaussian 
orbitals).  The  solid,  dotted,  and  dash-dotted  lines  in  Fig.  1  present  the 
calculated  energy  levels  of  s-,  p-,  and  d-states  of  hydrogen  atom,  respectively. 
As  the  numbers  of  the  Gaussian  functions  in  the  basis  sets  increase  from  5 
to  15,  the  calculated  energy  levels  approach  their  converged  values.  The 
calculated  first  three  energy  levels  of  hydrogen  atom  converged  to  -1.0,  - 
0.250,  and  -0.111  Ry,  when  we  utilized  15  -  60  Gaussian  functions  in  the 
basis  sets.  These  calculated  energy  levels  of  hydrogen  atom  agree  very  well 
with  that  of  the  exact  solution  for  the  principal  quantum  number  n  =  1,  2, 
and  3,  respectively.  The  lowest  energy  level  in  Fig.  1  belongs  to  the  Is  state 
of  H  atom.  The  second  lowest  level  includes  2s  and  2p  states.  The  third  one 
includes  3s,  3p,  and  3d  states. 

When  the  numbers  of  the  utilized  Gaussian  functions  in  the  basis  sets 
are  larger  than  about  70,  the  calculated  electron  energy  levels  of  the  3p  and 
3d  states  of  hydrogen  atom  dropped  to  lower  values  than  those  of  the  exact 
solution.  The  features  illustrate  the  instabilities  of  the  variational  calcu- 
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lations  of  the  Ritz-process.  Fig.  1  also  shows  that  the  numerical  property 
of  the  instability  in  the  computations  of  the  3d  states  of  hydrogen  atom  is 
different  from  that  of  the  3p  state.  The  calculated  energy  level  of  the  3d 
state  of  hydrogen  atom  dropped  rapidly  when  the  numbers  of  the  utilized 
Gaussian  functions  in  the  basis  sets  are  larger  than  about  70.  The  calcu¬ 
lated  energy  level  of  the  3p  state  decreases  much  slowly,  as  compared  with 
that  of  the  3d  state  in  the  same  computations.  The  calculated  energy  level 
of  3s  state  remained  the  same  as  that  of  the  exact  solution  in  these  test 
computations. 

3.2  Electronic  Properties  of  ZnSe 

Zinc  selenide  (ZnSe)  is  a  member  of  the  II- VI  family.  Its  crystal  structure 
belongs  to  the  space  group  of  Tj.  In  these  calculations,  we  used  a  lattice 
constant  of  a  =  5.65A. [33 ,  34,  39]  The  formation  of  the  Zn-Se  bond  has  a 
partial  ionic  character  in  addition  to  its  covalent  bond.  We  first  calculated 
the  atomic  wave  functions  of  Zn°  and  Se°.  We  utilized  these  atomic  wave 
functions  as  the  basis  set  to  perform  the  ab-initio  calculation  of  the  electronic 
structure  of  ZnSe.  The  charge  transfer  in  ZnSe  was  estimated  from  the 
calculated  electronic  wave  functions.  A  Zn  atom  loses  about  1.1  electrons  to 
a  Se  atom.  We  then  calculated  the  atomic  wave  functions  of  the  ionic  states 
of  Zn1+  and  Sel~ ,  utilizing  self-consistent  ab-initio  atomic  calculations.  The 
radial  parts  of  the  atomic  wave  functions  were  expanded  in  terms  of  Gaussian 
functions.  A  set  of  even-tempered  Gaussian  exponentials  was  employed  with 
a  minimum  of  0.15  and  a  maximum  of  1.770  x  105.  We  used  18  Gaussian 
functions  for  the  s  and  p  states  and  15  for  the  d  states  of  Zn.  We  utilized  19 
Gaussian  functions  for  the  s  and  p  states  and  17  for  the  d  states  of  Se.  In  the 
self-consistent  ab-initio  calculations  for  the  electronic  structure  of  the  solid 
state  ZnSe,  we  included  a  mesh  of  28  k-points,  with  proper  weights  in  the 
irreducible  Brillouin  zone.  We  took  the  Is,  2s,  and  2p  states  of  Zn  and  Se  as 
the  core  states  in  the  frozen  core  approximation,  and  allowed  all  other  states 
to  relax  in  the  self-consistent  ab-initio  calculations.  The  computational  error 
for  the  charge  density  was  about  0.0023  for  44  electrons.  The  self-consistent 
potentials  converged  to  a  difference  around  10-5  after  about  60  iterations. 
The  total  number  of  iterations  varies  with  the  input  potentials. 

We  performed  five  distinct  calculations  in  order  to  determine  the  op¬ 
timal  basis  set  for  the  calculation  of  the  electronic  structure  of  ZnSe,  uti¬ 
lizing  the  procedure  as  discussed  in  the  above  section  and  in  the  previous 
publications. [30,  31]  Table  1  illustrates  the  atomic  orbitals  of  the  basis  sets 
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of  the  five  self-consistent  calculations.  We  first  carried  out  completely  self- 
consistent  calculations  for  ZnSe  using  a  basis  set  that  included  Zn(ls2s3s4s 
2p3p  3d)  and  Se(ls2s3s4s  2p3p4p  3d)  orbitals.  We  then  repeated  the  self- 
consistent  calculation  using  the  above  basis  set  as  augmented  by  the  orbitals 
describing  the  next  lowest-lying  energy  level  of  Zinc.  Hence,  Zn(4p)  orbitals 
were  added  to  basis  set  I.  We  then  plotted  the  energy  bands  from  these 
two  distinct  calculations,  i.e.,  calculation  I  and  II,  and  compared  them  nu¬ 
merically  and  graphically.  Differences  were  obvious.  We  then  performed 
calculation  III  and  compared  its  results  to  those  of  calculation  II.  Up  to  five 
self-consistent  calculations  were  performed  and  the  results  of  each  calcula¬ 
tion  were  compared  to  those  of  the  previous  one  as  explained  above. 

In  Fig.  2,  we  presnt  the  calculated  electron  energy  bands  of  ZnSe  from 
calculation  III  (solid  lines)  and  calculation  IV  (dashed  lines).  Here  F  = 
(0, 0, 0);  L  =  (1, 1,  l)7r/a;  X  =  (0, 1, 0)2tt/o;  K  =  (3/4, 3/4, 0)27r/a.[35]  The 
zero  of  the  energy  was  set  at  the  top  of  the  valence  bands.  The  calculated 
energies  of  the  occupied  valence  bands  and  the  unoccupied  conduction  bands 
converged  to  an  average  difference  of  several  meV.  The  fourth  calculation 
(the  dashed  line)  of  Fig.  2  gives  sufficiently  converged  electron  energy  bands 
of  ZnSe  with  respect  to  the  size  of  the  basis  set.  As  we  present  in  the  follow¬ 
ing  sections,  the  calculated  electronic  properties  of  ZnSe  from  calculation 
IV  agree  well  with  experimental  data. 

However,  as  we  added  more  atomic  orbitals  in  the  calculations,  the  cal¬ 
culated  electron  energy  bands  of  ZnSe  drifted  away  from  the  converged  re¬ 
sults  of  calculation  IV  around  the  F-point.  One  may  observe  a  computa¬ 
tional  behavior  that  is  similar  to  that  of  the  calculations  for  hydrogen  atom. 
Fig.  3  shows  the  comparison  of  results  from  calculation  IV  (solid  lines)  and 
from  calculation  V  (dashed  lines).  The  calculated  occupied  energy  levels 
converged  with  respect  to  the  size  of  the  basis  set.  However,  some  of  the 
calculated  unoccupied  energy  levels  become  lower  than  those  of  the  calcu¬ 
lation  IV,  as  more  atomic  orbitals  are  added  in  the  calculations.  Different 
from  the  case  of  hydrogen  atom,  there  was  no  an  exact  solution  for  the 
electronic  structure  of  ZnSe  as  a  reference.  Various  physical  interpretations 
were  devised  to  explain  the  numerical  performance  of  the  Ritz-process  in 
the  calculations.  Some  of  them  were  believed  to  contribute  in  various  as¬ 
pects  to  the  problem.  In  fact,  the  problem  was  previously  claimed  as  the 
failure  of  density  functional  theory  or  the  failure  of  local  density  functional 
calculations  in  many  occasions.  Part  of  the  support  for  this  claim  was  due 
to  the  fact  that  the  calculated  conduction  band  energies  (or  the  band-gap 
energies)  of  semiconductors  utilizing  various  computational  methods  often 
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were  much  lower  than  that  of  experimental  results. 

In  spite  of  the  fact  that  the  mathematical  properties  of  different  compu¬ 
tational  methods  are  different,  it  can  be  proved  that  the  basis  sets  used  in 
those  of  the  previous  ab-initio  calculations  did  not  satisfy  the  requirement 
of  the  strongly  minimal  systems.  Those  include  all-electron  LAPW  cal¬ 
culations,  LMTO  calculations,  pseudopotentials  orthogonalized  plane  wave 
methods,  and  ab-initio  Gaussian  orbitals  or  LCAO  calculations.  Briggs  et 
al  used  multi-grid  accelerations  and  a  real-space  grid  as  a  basis  in  their  cal¬ 
culations  of  the  electronic  structure  of  elongated  diamond  cell,  an  isolated 
C60  molecule,  and  a  32-atom  cell  of  GaN.[36]  As  Briggs  et  al  discussed  in 
their  article,  the  real-space  grid  implementation  can  also  introduce  a  spu¬ 
rious  dependence  of  the  Kohn-Sham  eigen-values,  the  total  energy,  and  the 
ionic  forces  on  the  positions  of  the  ions  with  respect  to  the  real-space  grid. 

Nevertheless,  the  mathematical  and  physical  properties  of  the  ab-initio 
electronic  structure  calculations  should  be  carefully  examined.  As  we  pre¬ 
sented  in  the  previous  sections,  the  utilization  of  the  non-strongly  minimal 
systems  of  the  basis  sets  also  play  a  role  for  the  problem.  In  principle,  an 
exact  solution  cannot  be  reached  as  the  size  n  of  the  basis  set  goes  to  infin¬ 
ity.  The  Ritz-coefRcients  may  scatter  significantly  from  one  approximation 
to  the  next  as  n  becomes  large.  However,  the  calculation  of  the  occupied 
states  has  to  satisfy  the  restriction  of  Eq.  (3)  which  may  restrict  the  scatter¬ 
ing  of  the  Ritz-coefficients  of  these  states.  The  calculation  of  the  unoccupied 
states  do  not  have  the  restriction  and  the  problem  of  the  scattering  of  the 
Ritz-coefficients  is  particularly  severe  when  the  sizes  of  the  basis  sets  be¬ 
come  relatively  large.  We  have  proposed  that  an  optimum  basis  set  may  be 
needed  so  that  (a)  the  calculated  electron  density  is  converged,  and  (b)  the 
significant  scattering  of  the  Ritz-coefficients  is  avoided.  In  the  calculation 
of  the  electronic  structure  of  ZnSe,  we  selected  the  basis  set  IV  as  the  opti¬ 
mal  one  that  lead  to  the  converged  electron  density  and  related  properties. 
The  calculated  electron  energy  bands  are  illustrated  by  the  dashed  line  of 
Fig.  2  or  by  the  solid  lines  of  Fig.  3.  One  may  observe  that  the  calculated 
energy  bands  of  the  larger  basis  set  (Calculation  V  in  Fig.  3)  are  usually 
lower  or  unchanged  as  compared  to  those  of  the  smaller  basis  set  (Calcu¬ 
lation  IV  in  Fig.  3).  This  property  can  be  understood  from  the  Rayleigh 
theorem.[41,  42,  30,  31] 

As  shown  in  Fig.  2  and  Fig.  3,  ZnSe  is  a  direct  band  gap  semiconductor 
with  the  smallest  gap  energy  of  2.6  eV  at  the  center  of  the  Brillouin  zone 
(r).  Our  calculated  band-gap  energy  is  in  good  agreement  with  experimental 
measurements  of  about  2.8  eV.[34,  37,  38,  39]  We  present  the  total  density 
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of  states  (DOS)  of  ZnSe  in  Fig.  4,  as  obtained  from  calculation  IV.  The 
total  DOS  curve  in  Fig.  4,  particularly  the  inset,  indicates  that  the  small 
values  of  the  density  of  states  from  2.6  eV  to  2.9  eV  strongly  suggest  that 
measurements  or  related  analysis  that  are  not  very  sensitive  may  not  detect 
the  smallest  energy  gap. 

Our  calculated  effective  mass  of  n-type  carriers,  m* ,  near  the  F-point  is 
0.15  mo,  where  mo  is  the  free-electron  mass.  This  result  is  in  good  agreement 
with  the  experimental  value  of  0.16mo.[40,  39]  The  calculated  effective  mass 
of  n-type  carriers,  away  from  the  F-point,  is  between  0.15mo  and  0.17mo. 

We  also  performed  the  self-consistent  ab-initio  calculations  for  the  elec¬ 
tronic  structure  of  ZnSe,  utilizing  local  density  approximation  for  the  exchange- 
correlation  potential.  We  compared  the  electron  energy  bands  from  the  LDA 
and  GGA  calculations.  The  average  difference  between  the  resulted  bands 
of  the  two  calculations  is  of  the  order  of  1  mRy. 

3.3  Electronic  Structure  of  Ge 

We  performed  the  self-consistent  ab-initio  calculations  for  the  electronic 
structure  of  crystalline  Ge  in  the  diamond  structure.  We  used  a  zero  temper¬ 
ature  lattice  constant  of  a  =  5.63A  in  these  calculations. [39]  We  tested  the 
convergency  of  the  electronic  structure  calculation  with  respect  to  the  size 
of  the  basis  sets  of  atomic  orbitals.  We  calculated  the  atomic  wave  function 
using  a  separate  computer  program  that  employed  ab-initio  density  func¬ 
tional  computations  for  atoms.  We  expanded  the  atomic  wave  functions  as 
a  linear  combination  of  Gaussian  orbitals  (LOGO)  in  real  space  as  discussed 
in  the  previous  section. 

In  the  calculation  of  the  electronic  structure  of  Ge,  we  employed  non¬ 
local  density  functional  potentials  in  the  generalized  gradient  approximation 
that  was  developed  by  Perdew  and  Wang.  [22]  We  first  carried  out  com¬ 
pletely  self-consistent  calculations  for  crystalline  Ge  using  a  basis  set  that 
included  atomic  orbitals  of  Ge(ls2s3s4s  2p3p4p  3d).  We  then  repeated  the 
self-consistent  calculation  using  an  augmented  basis  set  that  also  included 
the  orbital  of  Ge(4d°).  In  the  self-consistent  calculations  of  the  electronic 
structure  of  solid  Ge,  we  considered  Ge(ls2s  2p)  as  the  core  states.  All  other 
states  were  treated  as  valence  states  and  were  allowed  to  relax  in  the  self- 
consistent  calculations.  We  then  plotted  the  resulting  energy  bands  from 
these  two  distinct  self-consistent  calculations.  We  observed  that  the  occu¬ 
pied  and  unoccupied  bands  from  the  two  calculations  differ  considerably. 

The  next  step  was  to  repeat  our  procedure  for  the  third,  fourth,  and 
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fifth  time.  In  these  calculations,  more  atomic  orbitals  were  added  one  at 
a  time.  The  results  of  the  fourth  calculation  (solid  lines)  and  of  the  fifth 
calculation  (dashed  lines)  are  shown  in  Fig.  5,  where  the  top  of  the  valence 
bands  is  set  as  the  zero  of  the  energy.  The  fourth  calculation  included  the 
atomic  orbitals  of  Ge(ls2s3s4s5s°  2p3p4p5p°  3d4d°).  The  fifth  calculation 
included  an  additional  atomic  orbital  of  Ge(5d°).  The  calculated  electron 
energy  bands  from  the  fourth  calculation  do  not  have  a  noticeable  difference 
from  that  of  the  fifth  calculation  for  a  large  energy  range  from  -13  eV  to 
about  8  eV.  The  average  difference  between  the  occupied  energy  bands  from 
the  two  calculations  in  Fig.  5  is  of  the  order  of  several  meV.  The  fourth 
calculation  (the  solid  lines  of  Fig.  5)  leads  to  sufficiently  converged  electron 
energy  bands  of  crystalline  Ge  with  respect  to  the  size  of  the  basis  set. 

We  present  the  calculated  energies  (eV)  of  the  electronic  states  of  Ge 
at  the  selected  high  symmetry  points  in  Table  2,  along  with  some  experi¬ 
mentally  measured  results.  ”Deg”  in  Table  2  refers  to  the  degeneracy  of  the 
electronic  states  at  that  k-point.  Our  calculated  results  for  the  electronic 
energies  agree  very  well  with  experimental  values.  The  calculated,  highest 
energy  state  of  the  valence  bands  is  at  the  T-point.  The  lowest  energy  state 
of  the  conduction  band  is  at  the  L-point.  The  calculated,  indirect  band  gap 
is  0.62  eV,  which  agree  well  with  experimental  results  of  0.66  -  0.74  eV.[39] 
In  these  calculations,  we  did  not  include  the  relativistic  effects  such  as  the 
spin-orbital  interaction  which  will  split  the  degenerate  bands  at  the  high 
symmetry  points. 

We  listed  the  calculated  effective  masses  of  n-type  carriers  (electrons) 
in  Table  3.  Although  the  spin-orbital  interaction  only  splits  the  electronic 
energy  levels  by  a  relatively  small  amount  (of  the  order  of  about  20  meV), 
the  curvature  (dispersion)  of  the  electron  energy  bands  can  change  more 
noticeably.  The  calculated  effective  masses  of  p-type  carriers  (holes),  with¬ 
out  the  consideration  of  the  spin-orbital  interaction  in  the  calculations,  can 
have  a  relatively  large  error  and  are  not  included  in  the  Table  3.  The  lowest 
conduction  band  is  not  degenerate  and  the  spin-orbital  interaction  may  not 
affect  substantially  the  effective  masses  of  n-type  carriers  in  this  band.  Our 
calculated  effective  masses  of  n-type  carriers  agree  well  with  experimental 
results.  [39] 

We  also  compared  the  electron  energy  bands  of  Ge  from  the  LDA  and 
GGA  calculations.  The  average  difference  between  the  energy  bands,  result¬ 
ing  from  the  two  potentials,  was  approximately  1  mRy. 


12 


3.4  Electronic  Structure  of  GaAs 

We  also  applied  the  same  calculation  procedure  as  discussed  in  the  previ¬ 
ous  sections  to  the  study  of  the  electronic  structure  of  GaAs  in  the  zinc 
blende  structure.  We  used  the  experimental  lattice  constant  of  a  =  5.654  ‘ 
in  the  calculations.  [43,  44]  We  calculated  the  atomic  wave  functions  of  Ga 
and  As  self-consistently.  We  evaluated  the  charge  transfer  in  GaAs  from  the 
self-consistent  ab-initio  calculations  of  the  electronic  structure.  Ga  atoms 
in  GaAs  lose  about  0.3  to  0.4  electrons  per  atom  to  As  atoms.  In  the  self- 
consistent  ab-initio  calculations,  we  employed  non-local  density  functional 
potentials  in  the  generalized  gradient  approximation.  We  used  the  Is,  2s, 
and  2p  states  of  Ga  and  As  as  the  core  states  in  the  frozen  core  approxi¬ 
mation,  and  allowed  all  other  states  to  relax  in  the  self-consistent  ab-initio 
calculations.  In  Fig.  6,  we  present  the  calculated  electron  energy  bands 
of  GaAs  along  some  high  symmetry  directions  in  the  Brillouin  zone.  Here 
T  =  (0,0,0);  L  =  (1, 1,  1)tt/o;  X  =  (0,  l,0)27r/a;  K  =  (3/4, 3/4, 0)27r/a. 
The  zero  of  the  energy  was  set  at  the  top  of  the  valence  bands.  In  Fig.  6, 
the  solid  lines  represent  the  calculated  results  using  the  atomic  orbitals  of 
Ga(ls2s3s4s  2p3p4p  3d4d)  and  As(ls2s3s4s  2p3p4p  3d4d),  where  Ga(4d) 
and  As(4d)  are  the  unoccupied  atomic  orbitals.  The  dashed  lines  represent 
the  calculated  bands  for  the  basis  set  that  has  an  additional  atomic  orbital 
of  As(5s).  The  calculated  electron  energy  bands  of  the  occupied  states  of 
GaAs  from  the  two  different  calculations  agree  very  well.  The  average  dif¬ 
ference  of  the  energies  of  the  occupied  states  from  these  two  calculations  is 
about  several  meV. 

GaAs  has  a  direct  band-gap  at  the  F-point  as  shown  in  Fig.  6.  The  lowest 
energy  state  of  the  conduction  band  is  at  the  F-point.  The  lowest  energy  of 
the  conduction  band  at  the  L-point  is  higher  than  that  at  the  F-point  by  0.2 
eV.  This  result  agrees  well  with  the  experimental  value  of  0.2  to  0.3  eV.[44] 
The  theoretical  band-gap  energy  of  GaAs  is  1.24  eV  from  this  calculation. 
The  reported  experimental  values  of  the  band-gap  energy  of  GaAs  is  1.4  to 

1.5  eV.[43,  44,  45]  As  we  discussed  in  the  previous  sections,  we  present  below 
the  calculated  electron  density  of  states  of  GaAs  in  Fig.  7.  There  is  a  tail 
structure  near  the  conduction  band  edge.  The  small  values  of  the  density  of 
states  from  1.24  eV  to  1.7  eV  strongly  suggest  that  measurements  or  related 
analyses  that  are  not  very  sensitive  may  not  detect  the  smallest  energy  gap. 

We  also  calculated  the  effective  mass  m*  of  the  electrons  (n-type  carriers) 
in  the  lowest  conduction  band  near  the  T-point.  The  effective  mass  m* 
is  calculated  from  the  curvature  of  the  band  dispersions.  The  calculated 
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effective  mass  of  the  electrons  do  not  show  much  difference  in  the  three 
directions  that  include  L  -t  T,  X  -t  f,  and  K  ->  T.  They  are  0.067, 
0.068,  and  0.068  mo,  respectively.  The  average  value  of  the  effective  mass 
of  the  electrons  is  m*  =  0.068mo,  where  mo  is  the  free  electron  mass.  The 
experimental  value  of  the  effective  mass  m*  of  the  electrons  is  about  0.065  to 
0.069  mo- [44]  Our  calculated  results  agree  well  with  the  experimental  data. 

We  compared  the  electron  energy  bands  of  GaAs  from  the  LDA  and 
GGA  calculations.  The  average  difference  between  the  energy  bands  of  the 
LDA  and  GGA  calculations  is  of  the  order  of  1  mRy. 

4  Conclusion 

In  this  article,  we  reexamined  some  of  the  mathematical  and  physical  prop¬ 
erties  of  the  ab-initio  density  functional  calculations  for  the  electronic  struc¬ 
ture  of  ZnSe,  Ge,  and  GaAs.  We  also  performed  test  computations  for  the 
electronic  structure  of  hydrogen  atom.  The  utilizations  of  non-strongly  min¬ 
imal  systems  in  the  self-consistent  ab-initio  LCAO  (or  LOGO)  calculations 
could  lead  to  the  non-uniformity  in  approaching  the  solutions  in  the  Ritz- 
process.  We  are  hot"  aware  a  previous  research  article  that  fully  addressed 
these  problems.  We  have  utilized  a  global  self-consistent  approach  to  de¬ 
termine  the  optimal  basis  sets  for  the  calculation  of  the  electronic  structure 
of  ZnSe,  Ge,  and  GaAs.  Our  calculated  results  of  the  electronic  properties 
agree  well  with  experimental  data. 
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Table  1:  The  atomic  orbitals  in  the  basis  sets  of  the  five  distinct  calculations 
for  the  electronic  structure  of  ZnSe. 


Basis 

Zn  valence 

Se  valence 

Bp 

ls2s  2p 

ls2s  2p 

3s4s  3p  3d 

3s4s  3p4p  3d 

ls2s  2p 

ls2s  2p 

3s4s  3p4p  3d 

3s4s  3p4p  3d 

ls2s  2p 

ls2s  2p 

3s4s  3p4p  3d 

3s4s  3p4p  3d4d  j 

IV 

ls2s  2p 

ls2s  2p 

3s4s  3p4p  3d4d 

3s4s  3p4p  3d4d  ; 

V 

ls2s  2p 

ls2s  2p 

3s4s  3p4p  3d4d 

3s4s5s  3p4p  3d4d 
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Table  2:  Calculated  electronic  energies  (eV)  of  Ge  at  the  selected  high  sym¬ 
metry  points.  The  experimental  data  are  from  Reference  [39,  44]. 


Ge 

Calculation 

Deg 

Measurements 

ri* 

-12.59 

1 

-12.6,  -12.9T0.2 

0.0 

3 

0.0 

T2c 

0.83 

1 

0.89 

Tisc 

2.81 

3 

3.01 

L'2v 

-10.41 

-10.6T0.5 

Liv 

-7.44 

■ 

-7.7±0.2 

-1.39 

2 

-1.4±0.3 

Lie 

0.62 

1 

0.74 

Lzc 

4.02 

2 

4.3±0.2,  4.2±0.1 

l'2c 

8.29 

1 

7.8±0.6,  7.9T0.1 

Xlv 

-8.64 

2 

-9.3±0.2 

Xiv 

-3.00 

2 

-3.15T0.2,  -3.5T0.2 

Xlc 

1.04 

2 

1.3±0.2 

Xzc 

9.98 

2 

Ea 

0.62 

0.66,  0.74 

Table  3:  The  effective  masses  (in  mo)  of  the  n-type  carriers  at  the  lowest 
conduction  band  of  Ge  in  the  diamond  structure. 


calculation 

measurements 

0.0807(39] 

| 

1.43 

1.57(39] 

IMWi! 

0.043 

0.038(39] 
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Figure  1:  The  calculated  electronic  energy  levels  of  hydrogen  atom  versus 
the  total  numbers  of  the  utilized  Gaussian  functions.  The  solid,  dotted,  and 
dash-dotted  lines  present  the  calculated  energy  levels  of  s-,  p-,  and  d-states 
of  hydrogen  atom,  respectively. 
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Figure  2:  The  calculated  electronic  energy  bands  of  ZnSe  along  the  high 
symmetry  directions.  The  solid  and  dashed  lines  represent  the  results  of 
calculations  III  and  IV,  respectively. 
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Figure  3:  The  calculated  electronic  energy  bands  of  ZnSe  from  Calculation 
IV  (solid  lines)  and  Calculation  V  (dashed  lines) . 
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Figure  4:  The  calculated  density  of  states  of  ZnSe.  The  inset  shows  the  tail 
structure  of  the  density  of  states  of  ZnSe  near  the  bottom  of  the  conduction 
bands. 


Figure  5:  The  calculated  electronic  energy  bands  of  Ge  along  the  high  sym¬ 
metry  directions.  The  solid  and  dashed  lines  represent  the  results  of  the 
fourth  and  fifth  calculations,  respectively. 
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Figure  6:  The  calculated  electron  energy  bands  of  GaAs. 
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Figure  7:  The  calculated  density  of  states  of  GaAs.  The  inset  shows  the  tail 
structure  of  the  density  of  states  of  GaAs  near  the  bottom  of  the  conduction 
bands. 
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